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^ Abstract 

A renormalization group (RG) improvement of the Einstein-Hilbert action is 
performed which promotes Newton's constant and the cosmological constant to 
scalar functions on spacetime. They arise from solutions of an exact RG equation 
by means of a "cutoff identification" which associates RG scales to the points of 
spacetime. The resulting modified Einstein equations for spherically symmetric, 
static spacetimes are derived and analyzed in detail. The modifications of the 
Newtonian limit due to the RG evolution are obtained for the general case. As an 
application, the viability of a scenario is investigated where strong quantum effects in 
the infrared cause Newton's constant to grow at large (astrophysical) distances. For 
two specific RG trajectories exact vacuum spacetimes modifying the Schwarzschild 
metric are obtained by means of a solution-generating Weyl transformation. Their 
possible relevance to the problem of the observed approximately flat galaxy rotation 
curves is discussed. It is found that a power law running of Newton's constant with 
a small exponent of the order 10~ 6 would account for their non-Keplerian behavior 
without having to postulate the presence of any dark matter in the galactic halo. 
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1 Introduction 



Renormalization group improvement is a standard tool of elementary particle physics. 
Typically it is used for including higher order quantum corrections into scattering am- 
plitudes or cross-sections by replacing the classical or "bare" parameters on which they 
depend with scale dependent "running" coupling constants and masses obtained by solv- 
ing an appropriate renormalization group equation. A nontrivial step in this procedure 
is identifying the physical scale (momentum transfer, etc.) at which the running cou- 
plings are to be evaluated. While originally employed within perturbation theory, this 
idea can also be used in the context of the various "exact renormalization group equa- 
tions" [1] which implement the idea of the Wilsonian renormalization group in quantum 
field theory. A particularly convenient and powerful formulation is the "effective average 
action" [2, 3] and its associated flow equation. One of its advantages is that it can be 
applied to gauge theories [4] and also to gravity [5] . 

In the latter case the effective average action Tk [g^u] is a free energy functional 
which depends on the metric and a momentum scale k with the interpretation of a variable 
infrared (IR) cutoff. The action T k is similar to the ordinary effective action T which it 
approaches for k — > 0. The main difference is that the path integral defining Y k extends 
only over quantum fluctuations with covariant momenta p 2 > k 2 . As a result, T k describes 
the dynamics of metrics averaged over spacetime volumes of linear dimension k~ x . The 
functional T k [g^] gives rise to an effective field theory valid near the scale k. Hence, 
when evaluated at tree level, T k correctly describes all quantum gravitational phenomena, 
including all loop effects, provided the typical momentum scales involved are all of order 
k. Considered a function of k, T k describes a renormalization group (RG) trajectory in 
the space of all action functionals. The trajectory can be obtained by solving an exact 
functional RG equation. In practice one has to resort to approximations. Nonperturbative 
approximate solutions can be obtained by truncating the space of action functionals, 
i. e. by projecting the RG flow onto a finite-dimensional subspace which encapsulates 
the essential physics. Recently a great deal of work went into the investigation of such 
truncated RG flows of quantum gravity [5-13]. Most of these studies employ the "Einstein- 
Hilbert truncation" which approximates T k by a linear combination of the monomials 
j^fg R and Jyfg. Their prefactors contain the running Newton constant G(k) and the 
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running cosmological constant A(k). Their ^-dependence is governed by a system of two 
coupled ordinary differential equations. 

Let us assume we found some solution to this system, a RG trajectory (G(k), A(k)Y 
How can we take advantage of this information? Sometimes it is possible to "RG improve" 
the predictions of classical General Relativity by replacing the pertinent classical (bare) 
parameters with the running ones, G(k) and A(k). In a spacetime formulation, one would 
try to find a "cutoff identification" of the form [14] 

k = k(x) (1.1) 

which converts the scale dependence of G and A to a position dependence: 

G{x) = G(k = k(x)), A(x) = A(k = k(x)). (1.2) 

Clearly it is not always possible to associate cutoff values to spacetime points and to 
find a map x i— > k(x). As we discussed in ref. [15], a detailed discussion of the possible 
physical cutoff mechanisms which might be operative in a given system is required here. 
However, typically when the problem has a high degree of symmetry, maps of the kind 
k = k(x) actually do exist and can be guessed almost unambiguously [15]. For instance, 
in Robertson- Walker cosmology the condition of homogeneity and isotropy leads to strong 
constraints on the possible form of k(x): it can depend on the cosmological time t only. 
In fact, the most natural choice is [16-19] 

kit) = £/*, (1.3) 

with £ a constant of order unity. The motivation for (jl.3|) is obvious: By the time the age 
of the Universe is t, no quantum fluctuation with a frequency smaller than 1/t can have 
played any role yet. Hence the integrating-out of field modes should be stopped at the 
IR cutoff k « 1/t. (See also [20,21].) 

Likewise, if the spacetime is spherically symmetric and static, k can depend on the 
radial coordinate r only. In ref. [22] the RG improvement of Schwarzschild black holes 
has been discussed and it was argued that the correct cutoff is 

k{r)=t/d{r), (1.4) 

where, outside the horizon, d(r) is the radial proper distance to the center of the black 
hole. 
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Knowing the position dependence of G{x) and A(x) the actual RG improvement 
can be performed in either of the following ways. (1) One picks a spacetime solving the 
classical Einstein equation and substitutes G — > G(x), A — > A(x) in this solution. (2) 
One performs the substitution one step earlier at the level of the field equations, and then 
solves the resulting "improved Einstein equation" 

= -A(x) 9ia , + 8vr G{x) T^. (1.5) 

(3) One performs the substitution still another step earlier, at the level of the action 
from which Einstein's equation is derived, then derives the field equations implied by the 
improved action and solves them. 

The three options in this list are ordered according to increasing predictive power 
and generality. The improved solution approach [22] can be used only when the quantum 
corrected metric is not very different from the classical one. The improved equation 
approach [17] does allow for large quantum corrections, but there are also situations where 
the improved Einstein equation ()1.5|) obviously does not capture the essential quantum 
effects. The prime example are vacuum spacetimes for which the cosmological constant is 
not an issue, small isolated systems such as black holes, for instance. While we expect on 
general grounds that the running of G is important for the structure of such spacetimes, 
the improved Einstein equation (jl.5|) cannot "know" about this running since it does not 
differ from the classical one, G^ v = 0. 

In ref. [17] the improved equation approach has been applied to cosmology and 
its relation to the improved solution technique was clarified. Ref. [22] contains the first 
application of the effective average action to black holes. It used the improved solution 
technique which is insufficient near the central singularity, for instance, where the quantum 
effects become strong. This was one of the motivations for developing the improved action 
approach in ref. [15]. In the present paper we are going to take a first step towards its 
application to black holes. 

The idea is to start from the classical Einstein-Hilbert action Seh = / d 4 x y/—g J2?eh 
with the Lagrangian =S?eh = (R — 2A)/ (167T G) and to promote, in this Lagrangian, G 
and A to scalar fields on spacetime. This leads to the modified Einstein-Hilbert (mEH) 
Lagrangian 

J? mEB =—^ 7 - S {R-2A(x)}. (1.6) 
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The resulting modified Einstein equations differ from the classical ones not only by the 
substitution G — > G{x), A — > A(x) but also by terms containing derivatives of G(x). 
Because of those terms, the improvement of the field equations and of the action functional 
are not equivalent a priori. In ref. [15] the relation between the two approaches was 
discussed in detail. The upshot of this analysis was that, in the improved action approach, 
the modified Einstein-Hilbert action S m EH = / d 4 x y/—g -5f m EH should be supplemented 
by an action Sg [g^ u , G, A] whose associated energy- momentum tensor 



describes the 4-momentum carried by the scalar fields G(x) and A(x). In the approach of 
improving field equations those scalars do not carry gravitating 4-momentum, but they do 
so in the approach based upon J2? m EH- The energy- momentum tensor Q^ v plays a double 
role: In physical situations where both approaches are applicable it renders their respective 
results equivalent, at least at a qualitative level, and at the same time it helps in making 
the modified field equations mathematically consistent. In fact, the modified Einstein 
equation has the form G^ = where contains the energy-momentum tensor of the 
matter system, T^ U) the scalars G(x) and A(x), and their derivatives. Since D^G^ = 
by Bianchi's identity, consistency requires that D^S^ = 0. In classical General Relativity 
this condition is satisfied if T^ v is conserved. In the RG improved setting additional 
conditions on A(x) and G(x) arise, henceforth referred to as "consistency conditions". 
As a rule, the consistency conditions become much less restrictive when an appropriate 
^-tensor is included. 

Before continuing, let us recall some properties of the RG trajectories of Quantum 
Einstein Gravity [8,9]. At the exact level, they consist of infinitely many generalized 
couplings G(k), A(k), ■ ■ ■ needed in order to specify a generic action function, including 
higher derivative and nonlocal invariants [12,23]. Rather than in terms of the dimensionful 
quantities G(k), A(k), ■ ■ ■ the flow equations are usually written in terms of their dimen- 
sionless counterparts g(k) = k 2 G(k), X(k) = A(k)/k 2 , Within the Einstein-Hilbert 
truncation, the RG equations for g{k) and X(k) were solved numerically in [8]. The RG 
flow on the g-A-plane is dominated by two fixed points (<?*,A*): a Gaussian fixed point 
at g* = A* = 0, and a non-Gaussian one with g* > and A* > 0. The high-energy, i.e. 
short- distance behavior of Quantum Einstein Gravity is governed by the non-Gaussian 
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fixed point: for k — > oo, all RG trajectories run into this fixed point. If it is present in 
the exact flow equation, it can by used to construct a fundamental, microscopic quantum 
theory of gravity by taking the limit of infinite UV cutoff along one of the trajectories 
running into the fixed point. The theory is said to be "nonperturbatively renormalizable" 
then [9,10,24,25]. At the non-Gaussian fixed point the dimensionless couplings assume 
constant values so that the dimensionful ones display a power law running: 

G(k) = gjk 2 , k(k) = Kk 2 . (1.8) 

Newton's constant G(k) is an asymptotically free coupling: it vanishes for k — > oo. 
Conversely, lowering k, one finds that G(k) increases monotonically all the way down 
from the ultraviolet (UV) to the IR. Interpreting 1/k as a distance scale, we see that 
G increases with distance: quantum gravity is anfecreening, a property it shares with 
QCD, but not with QED, for instance. For the RG trajectory ending at X(k = 0) = the 
scale dependence of G can be approximated by 

1 + wGr 

where a; is a constant. Here G = G(k = 0) is the renormalized Newton constant. If 
the behavior (jl.9j) applies, this is the value we measure in the laboratory because the 
difference between G at k = and at a typical laboratory scale, k — (1 meter) -1 say, is 
utterly negligible. According to (jl.9J) . G(k) is practically constant once k is smaller than 
the Planck mass mpi = 1/ \f~G . 

It is known [8] that the Einstein-Hilbert approximation performs better in the UV 
than in the IR. In particular it seems clear that for A; — ^ it should be supplemented by 
more complicated (nonlocal) invariants. One might wonder what their impact on the RG 
evolution in the IR could be. It is an exciting conjecture that gravity is subject to strong 
quantum corrections at large distances [26,27] which might lead to visible renormalization 
effects at large, i.e. astrophysical or cosmological distance scales [17,28,29]. From the 
theoretical point of view, this conjecture is supported by the observation that Quantum 
Einstein Gravity has many features in common with non-Abelian gauge theories [8,12]. 
In Yang-Mills theory, QCD, say, we know that the large-distance physics is governed 
by strong and highly nonperturbative quantum effects which are responsible for color 
confinement and chiral symmetry breaking, for instance. Because of the extreme algebraic 
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complexity of the corresponding flow equations it has not yet been possible to explore the 
IR of Quantum Einstein Gravity with realistic nonlocal truncations; see [12], however, for 
a first step in this direction. 

From the phenomenological point of view, the astrophysical and cosmological im- 
plications of running gravitational coupling constants were explored in [28, 29] in the 
framework of a higher-derivative theory of gravity; in particular, primordial nucleosyn- 
thesis, structure formation, and galaxy rotation curves were discussed. While the observed 
flatness of galaxy rotation curves v (r) [30] is usually considered evidence for the existence 
of dark matter, it was argued that macroscopic quantum effects could mimic the presence 
of dark matter. 

Similarly, the phenomenological theory of the MOdified Newtonian Dynamics called 
"MOND" [31] is designed so as to explain the galaxy rotation curves without the need 
for dark matter. For a point mass M, say, its modified (nonrelativistic) field equation 
reproduces the familiar — GM/r-potential at small distances; at large distances it yields 
instead 



which implies a flat rotation curve. Here Oq ~ 10 _8 cmsec -2 is a universal constant. If 
it exists, the In r-term in fll.lOj) overwhelms the classical 1 / r-potential at large distances. 
In quantum gravity the perturbative corrections to Newton's potential [32-35] vanish for 
r — > oo, so a lnr-behavior would require strong nonperturbative effects. 

A fully relativistic theory which also tries to explain the rotation curves as due to a 
non- Newtonian force is the Weyl theory put forward by Mannheim [36]. Its Lagrangian 
is the square of the Weyl tensor, leading to a potential of the type 1/r + r. The resulting 
rotation curves are not flat, but seem to be in accord with the observations. 

For a review of further work on non-standard gravity at galactic scales and a com- 
prehensive list of references we refer to [37] . In this context we also mention the possibility 
that the flatness of the rotation curves is caused by the quintessence field [38] originally 
introduced for cosmological reasons [39,40]. 

As for strong (quantum or classical [41]) renormalization effects at cosmological 
scales, an infrared fixed point model was proposed in ref. [17]. Its basic assumption is 
that, for k — > 0, the RG trajectory of the Universe is attracted to a second fixed point 




(1.10) 
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(d™, ^ R ), implying G(k) = g^/k 2 and A(k) = A* R k 2 or, with k oc 1/t, G oc t 2 and 
A oc 1/t 2 at late times (t — ► oo). It was not yet possible to derive this IR fixed point from 
a flow equation. Its main motivation is the interesting phenomenology it entails [17, 18]. 
In particular it explains without any finetuning why in the present Universe the matter 
energy density equals (approximately) the vacuum energy density. At the same time it 
performs with respect to the standard tests (supernovae, compact radio sources, etc.) as 
well as the best-fit Friedmann model [18]. 

One of the central questions we are going to address in the present paper is whether 
the flat galaxy rotation curves can be explained as due to the onset of the cosmological 
fixed point running at galactic scales. We shall employ a model trajectory of the form 



which interpolates between a constant G at large k (small distances) and the l//c 2 -running 
at small k (large distances). Since the cosmological constant is not expected to be essential 
at the scale of galaxies we set A(fc) = right from the start. The ansatz (jl.llj) is assumed 
to be valid at distances larger than laboratory scales, G being the measured value of 
Newton's constant. This ansatz neglects the UV renormalization effects of (jl.9j) which 
are completely irrelevant for galaxies. 

Within the improved action approach, we shall find that the G(k) of (jl.llj) has 
indeed the right qualitative properties. The running of G does not only prevent the 
rotation curves from falling off at large distances, they show a significant growth even. 
While many of the observed rotation curves are indeed gently rising [30], it appears that 
the growth implied by (jl.llj) is too strong and the relevant velocities v(r) are too small. 

A second model trajectory which we are going to analyze assumes a power law 



over a wide range of scales. We shall see that if the observed flatness of the rotation 
curves is indeed a running-G effect, the ansatz (jl,12j) should come very close to the truth 
about the underlying trajectory. At sufficiently large distances it gives rise to perfectly 
flat rotation curves. 

The purpose of the present paper is twofold. On the theoretical side, we study the 
application of the improved action approach of ref. [15] to static isotropic (spherically 



G(k) =G + g^/k 2 



(1.11) 



G{k)ock- q , 0<g«l, 



(1.12) 
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symmetric) spacetimes for an arbitrary RG trajectory. We derive the corresponding mod- 
ified Einstein equations, analyze their general structure, and construct a cutoff which is 
selfconsistent in the following sense. Employing a cutoff identification k(r) oc l/d(r) one 
has to specify a metric by means of which the radial proper distance d(r) is to be com- 
puted. When the quantum effects are weak one can use a fixed background metric, but 
in general one would like to compute d(r) using the actual metric, i. e. the solution of the 
modified Einstein equation. Clearly the problem is that by the time one writes down the 
cutoff and the resulting modified field equation one does not know its solution yet. We 
shall solve this problem by a kind of "bootstrap" which allows us to compute the relation 
between k and the coordinate r in such a way that k(r) oc l/d g (r) where d g is the proper 
distance with respect to the solution g^ v . 

These results will be important for a later application of the improved action ap- 
proach to black holes, for instance, with a G(k) of the type (jl.9|) . 

In the present paper, on the applied side, we are going to use these methods in order 
to explore the impact a scale dependent Newton constant could have on the properties 
of (idealized) galaxies. We shall derive exact results for the model trajectories and 
()I.I2j) and allow for an arbitrary G(k) in the (modified) Newtonian approximation. 

With regard to the underlying RG trajectories the application of (jl.llj) and (Jl.f2j) 
to galaxies is certainly more speculative than the use of ()I.9|) in black hole physics since 
the latter G(k) has actually been computed within the Einstein-Hilbert truncation, while 
the former ones have a phenomenological motivation only. The reason why we neverthe- 
less start with the galaxies as a first application is that here the essential new physics 
can be understood in a quasi-Newtonian language. In fact, a large portion of the present 
paper (Section |7J) will be devoted to the modified Newtonian limit of the improved ac- 
tion approach. This part of the paper is completely general, no special symmetries or 
trajectories are assumed. 

The remaining sections of this paper are organized as follows. In Section |21 we review 
some properties of the improved action approach which will be needed later on. In Section 
El we discuss the modified Einstein equations in the static isotropic case, and in Section |U 
we use Weyl transformations in order to generate a large class of solutions and to solve the 
problem of the selfconsistent cutoff identification. Up to this point the discussion applies 
to any scale dependence G(k). Then, in Sections and El we construct and analyze the 
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static isotropic spacetimes resulting from the running G of the fixed point scenario, eq. 
((l.lip . and the power law (|1.12|) . respectively. In Section[7|we return to the general theory 
and explore its modified Newtonian limit. At the end of this section we describe a simple 
RG-based model of galactic haloes and demonstrate explicitly that, with an appropriate 
RG trajectory, no dark matter would be needed in order to explain the flatness of the 
observed rotation curves. Section |H] contains the conclusions. 



2 RG improved gravitational actions 

In this section we briefly review the improved action approach of ref. [15] and discuss 
some new aspects which will be important later on. We investigate gravity interacting 
with an arbitrary set of minimally coupled matter fields, A(x). Their dynamics is governed 
by the action Syt[g, A] which gives rise to the energy- momentum tensor 

T „ v = 2 SS M \g,A] 

Since Su is diffeomorphism invariant by assumption, T^ v is conserved when A satisfies 
its equation of motion, 5Su/SA = 0: 

DpT 1 ™ = 0. (2.2) 

The total action is taken to be 

Stot = S mE}l [9, G, A] + S u [g, A] + S e [g, G, A] (2.3) 
with the modified Einstein-Hilbert action 2 

As we mentioned already, the purpose of the term Sg[g,G,A] is to generate the energy- 
momentum tensor 9^ v of the scalar fields G(x) and A(x). A priori Sg is quite arbitrary 
but, as we shall see, the consistency condition will lead to strong constraints on Sg and 6^. 



2 Our curvature conventions are R a piiv = —d^ fip a + R^v = The metric signature is 

(- + ++). We abbreviate (DG) 2 = g^ v D^GD V G and D 2 G = g^" D p D v G. The labeling of the 
Schwarzschild coordinates is x° — t, x 1 — r, x 2 — 6, x 3 = ip. 
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The only general requirements are that S$ is independent of A, that it is diffeomorphism 
invariant if G and A are transformed as scalars, and that 9^ u vanishes for G,A = const 
(in which case the scalars carry no momentum). 

Varying S to t with respect to the metric we obtain the modified Einstein equation 

G^v = -A 9»u + 8nG (T^ + AT„ V + 9^) , (2.5) 

or, with the convenient definitions -d^ = 8tcG9^ u and At^ = 8nG AT^ U , 

Gfj, u = -A + 8nG + At^ + d^. (2.6) 

The tensors T^ v and 9^ v obtain as the functional derivatives of Su and S$, respectively, 
and AT^y originates from the x-dependence of Newton's constant in 5' m EH : 

At^ = G(x){D^D v -g^ v D 2 ) ' 



V ; (2.7) 
= -L {2 D„G D V G - G D,D V G - 9flv [2 (DG) 2 - G D 2 G] } . 

The modified Einstein equation is coupled to the equations of motion of the matter system, 

-Ja = 1a =0 - (2 - 8) 

There is, however, no corresponding field equation for the scalars G(x) and A(x). They 
are "external fields" defining the "background" for the dynamics of g^ v and A. 

Not every choice of G(x), A(x), and 9^ v leads to consistent field equations. Since the 
LHS of ()2.6|) has vanishing covariant divergence, consistency requires the same to hold 
true for its RHS. This implies the "consistency condition" 

Z^AV + - D V A + 8tt (D„G) T% = 0, (2.9) 

or, in terms of AT^ U and 9^ u , 

8nG D" [AT^ + 9^] + 8tt [T„ u + AT„ U + 9„ u ] D^G — D U A = 0. (2.10) 

In writing down ()2.9|) and ()2.10|) we used that D^T^ V = 0. After working out the covariant 
divergence of At^ u , eq. (|2.9|) becomes 

G- 1 D»G (AV - R pv ) + - D U A + 8tt (D^G) T% = 0. (2.11) 
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Let us fix a constant reference value of Newton's constant, G, and let us define the 



"total" energy-momentum tensor T*° according to 



8ttG T™ = -A giw + 8ttG (T^ + AT^ + V) • (2-12) 
Then the modified field equation assumes the form of the classical equation with A = 0: 



G 



(IV 



8ttGT, 



tot 

(IV ' 



(2-13) 



The consistency condition (|2.9jl is nothing but the continuity equation for T * ot which 
guarantees the integrability of (j2.13j) : 



rjMT"' tot 



0. 



(2.14) 



Occasionally it is helpful to employ the following alternative form of the consistency 
condition: 

3 D U G 



2 G 3 



GD 2 G-2 (DGY 



(IV 



G » v 



1 



(2.15) 



+ AnTD v G--D v (GA) = 0. 

Gr 



Here = fl^ — g^v ^cT/2 and T ee T ^. Eq. (|2.15|) follows from eq. (j2.11|) by making use 
of the modified Einstein equation. In fact, eq. (|2.15j) is equivalent to eq. (|2.9j) only "on- 
shell" , i. e. when g^ v satisfies its equation of motion. However, since the field equations 
anyhow have to be solved together with the consistency condition we may equally well 
use ()2.15|) rather than ()2.9|) as the independent equation. 

In order to get a first impression of the content of the "on-shell consistency condition" 
(|2.15|) let us assume that T^ u is traceless (T = 0), at least when A satisfies its equation 
of motion, and that the ^-dependence of G and A is such that (GA) = 0. In this case 
(|2.15|) boils down to 

3 D U G 
2~G r 



(IV 



G ^ 



GD 2 G-2 (DGf 



0. 



(2.16) 



In ref. [15] we showed that there exists a unique tensor = which solves f!2.16|) 
identically in G(x), i.e. without assuming any special properties of this function, and 
which vanishes for G = const. The corresponding 6^ D = $^ D /87rG reads 

3\ 1 
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BD 

(IV 



87TG 3 



D^G D V G - \g^ u (DGY 



(2.17) 
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The superscript "BD" stands for "Brans-Dicke" . In fact, the tensor (j2.17j) made its first 
appearance in standard Brans-Dicke theory [42] . For a singular choice of the Brans-Dicke 
parameter, u = —3/2, (j2.17|) is exactly the energy-momentum tensor of the BD-scalar 
</> = 1/G. In ref. [15] we described in detail the relation of our theory to Brans-Dicke 
theory. This discussion is not relevant for the purposes of the present paper. As the 
prototype of an admissible ^-tensor, 9^P will play an important role, though. It can be 
derived from the action 

S e BD [9,G] = J- |d 4 x v^^|^. (2.18) 

In ref. [15] we identified various classes of solutions which enjoy special properties. 
Here a "solution" is a set {g^ u , A, G, A, 6^ which satisfies the coupled system of equa- 
tions consisting of the modified Einstein equation, the matter equation of motion, and 
the consistency condition. The specification of a class involves (i) a choice of 8^, (ii) as- 
sumptions about the scalar fields G(x) and A(x), and (iii) assumptions about the matter 



sector. The defining properties of the classes considered are: 

Class I: 6^ = (2.19a) 

Class II: 9„ v = ^ D , A = 0, T = (2.19b) 

Class Ilia: 6^ = 9® D , GA = const ^ 0, T = (2.19c) 

Class Illb: 9„ v = 0™, A ^ 0, d„ (GA) ^ (2.19d) 



In Class II and Class Ilia the on-shell consistency condition (J2.15|) is satisfied by con- 
struction. In Class Illb it reduces to the simple requirement A.nG T d^G = <9 M (GA) ^ 0. 
The importance of the solutions in Class Ilia resides in the fact that, near a RG fixed 
point, G(k) oc 1/k 2 and A(A;) oc k 2 which, for every cutoff identification k = k(x), implies 
G(a;)A(a;) = const. 

In the cases where T = we assume that the tracelessness of T^ v is due to a 
Weyl-invariance of the matter action, i.e., if is the Weyl weight of the matter field, 
S , M[exp(2cr) g^ u , exp(— 2A A a) A~\ = Sm [g^, A] for any a(x). This condition insures 
that T = if 6S M /SA = 0. 

In ref. [15] we demonstrated that all solutions (g^, A) of Class II can be generated 
in the following very efficient way. Let us pick a solution (7^, .4.) of the conventional 
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Einstein equation with G = const = G, A = 0, coupled to the matter equation of motion: 



<V(7) =87tGT^(A7), (2.20a) 
5S M [7,„4]/<U= 0. (2.20b) 

Now we use the prescribed field G(x) in order to perform the following Weyl transforma- 
tion of the metric 7^ and the matter field A: 

^(x) = ^ 7 ^(x), (2.21a) 



A(x) 



G(x^- A 
G 



A 



A{x). (2.21b) 



It can be shown that the pair (g^, A) thus obtained is a Class II solution of the RG 
improved system, with a position dependent Newton constant G = G(x), and that all 
solutions of Class II can be obtained in this way. 

The proof makes essential use of the fact that, for T = and A = 0, the Brans-Dicke 
action S(P B can be absorbed into S m ER by Weyl-rescaling the metric according to ()2.21aj) . 
A Weyl invariant matter action is "blind" to this rescaling, explaining why the case T = 
is special. (See [15] for further details.) 

Generically a nonzero cosmological constant with an arbitrary x-dependence will 
spoil the simplicity of this method. An exception is the case when A varies according to 

A(x) = A G/G{x), A = const. (2.22) 

Hence AG = const and we are in Class Ilia. It turns out that all solutions (g^, A) of 
Class Ilia can be obtained as the Weyl transform according to ()2.21|) of some solution 
(l^u,A) of the classical equations 

G^( 7 ) = -A 7^ + SnGT^A^), (2.23a) 
6S M [l,A}/5A=0. (2.23b) 

For a proof and further details about the use of Weyl rescalings as a solution-generating 
transformation the reader is referred to ref. [15]. 
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Very often one would like to find metrics g^ u which enjoy special symmetry proper- 
ties. For a symmetry transformation generated by the Killing vector field the condition 
is that the Lie-derivative of g^ v with respect to vanishes: S^Kg^v = 0. If (j2.21a|) ap- 
plies, i.e. in Classes II and Ilia, this leads to the requirement G J^kI^v + I^u-^rG = 0. 
For this requirement to be met it is sufficient, but not necessary, that K is a Killing 
vector of 7^, too, ^kI^v = 0, and that the "external field" G(x) is chosen so as to 
respect the symmetry: S^kG = d^G = 0. In the present paper we shall always gen- 
erate symmetric g^s from T^'s and G's possessing the same symmetry. (In principle a 
more general situation is conceivable where the non-invariance of G(x) compensates for a 
non-invariance of 7^.) 

3 Static isotropic solutions 

3.1 Symmetry properties 

In this section we focus on solutions to the modified Einstein equations with special 
symmetry properties. We search for spacetimes which are static and isotropic, i. e. spher- 
ically symmetric. The corresponding g^ v admits 4 Killing vectors: one related to time 
translation, and 3 to spatial rotations. In Schwarzschild-type coordinates x M = (t, r, 9, if) 
the most general metric of this kind can be brought to the "standard form" [43] 

ds 2 = dx» dx v = -B(r) dt 2 + A(r) dr 2 + r 2 da 2 , (3.1a) 

where 

da 2 = d9 2 + sin 2 6 dip 2 (3.1b) 

denotes the line element on the unit 2-sphere. For the spacetimes discussed in the present 
paper the functions A(r) and B(r) are strictly positive and horizons are not an issue. 
Hence t and r have the interpretation of the time and the radial coordinate, respectively. 

In the following the "external fields" G(x) and A(x) are chosen in accord with the 
above symmetry requirements. We demand that d^G = d^A = for all four 
Killing vectors which implies that G and A can depend on the radial coordinate r 
only: 

G = G(r), A = A(r). (3.2) 
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Likewise the matter field configuration is assumed to respect all symmetries. 



3.2 Decomposition of the field equations 

We model the matter system by a perfect fluid in hydrostatic equilibrium. In the 
metric (j3.1|) its energy-momentum tensor T v = pg/ + (p + p) u^u" is diagonal, 

Tp" = diag[-p,p,p,p]. (3.3) 

The energy density p and the pressure p depend on r only. To fully characterize the fluid 
one has to specify an equation of state 

p = p(p). (3.4) 
For the metric (J3.1)) the continuity equation D^T u = amounts to a single condition: 

(3.5) 



d B'(r) 



dr v ' 2B(r) 

(Here and in the following a prime always denotes the derivative with respect to the 
argument.) 

Next let us evaluate the tensor AT^ for metrics of the "A5-form" (|3.1|) . Working 
out the covariant derivatives in (|2.7|) one finds that AT " is diagonal: 

AT/ = At;/ (8tt G) = diag[-Ap, A Pr , Ap ± , Ap ± ] . (3.6) 

Here Ap = — AT t 1 is an energy density which is induced by the position dependence of 
Newton's constant: 



Ap 



1 



8ttG 



r A V G 



+ 



A' ( G' 



2A 2 \ G 



+ 



A \ G 



G" 
A \~G 



Likewise there is a contribution to the radial pressure, 



Ap r = AT/ 
and to the tangential pressure: 
A P± = AT e e = AT J 



87TG 



B' 



2AB V G 



G' 



+ 



rA \~G 



(3.7a) 



(3.7b) 



8ttG 





'G'\ 




k g) 


2 




"I 





A' ( G' 



2A 2 \G 



1 / G" 



B' ( G' 



2AB \ G 



A \ G 



(3.7c) 
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As expected on symmetry grounds, the components AT e e and AT ^ are found to be equal. 
However, in general the radial and the tangential pressures are different so that AT " is 
not of the perfect fluid type. 

The precise form of the (9-tensor will be left open for the time being. It is enough 
to know that it has the structure 

0; = (SttG) = diag[-p e , pe, r , p ,±, p e ,x] . (3.8) 

It supplies additional contributions pg, pg ;r , and pg ; ± to the energy density and the radial 
and tangential pressure, respectively, which are functionally dependent on A(r) and B(r). 

A further source term appearing on the RHS of Einstein's equation is the one con- 
taining the cosmological constant. Setting T A ^ = — (A/8nG) g" = diag[— Pa, Pa, Pa, Pa] , 
the vacuum energy density and pressure due to the cosmological constant are, respectively, 
Pa = A/ (87rG) and p\ = —A/ (8nG). In this notation, the total energy- momentum tensor 
defined by eq. (j2.12|) reads T tot ^ = diag[— p tot ,p r tot ,p_J; ot ,pj|; ot ] with the entries 

p**= [G(r)/G] (p + PA + Ap + p e ), 

p r tot = [G(r)/G] (p + p A + Ap r +p e , r ), (3.9) 

pi ot =[G(r)/G] (p + PA + Ap ± +p e , ± ). 

For the metric (|3.1|) the only nonvanishing components of the Einstein tensor G" are 
G t \ G r r , and G e e = G*. As a result, the modified Einstein equation G£ = 8rrGT tot ^ 
decomposes into the following 3 equations: the tt-component, 

G t = - ^ - ^ = (p + PA + Ap + pg) , (3.10a) 



the rr-component, 



1 1 B 1 

G r = — A - - + —T5 = 87rG (P + Pa + APr + P9,r) , (3.10b) 
r 2 A r z r A B 



and the 99- or <p<p-component: 



2rA 2 2rAB AA 2 B AAB 2 2AB 
8nG (p + p A + Ap ± +p e , ± ). 



(3.10c) 
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In the same notation, the consistency condition assumes the form D^T 
Working out the covariant derivative leads to 



tot v 



0. 



% Pr t0t + ^(P + Pr t0t ) + ~ (P^ ~ P f* ) = 



(3.11) 



whence 



G' (p + Pa + Ap r + p e ,r) +G — (p A + Ap r + p 0ir ) 



B' 



(3.12) 




In deriving ()3.12j) we exploited that T^ u obeys the ordinary continuity equation (|3.5j) . 

As in the classical case, Einstein's equations are not independent of the conserva- 
tion law for the energy-momentum tensor. Applying the standard discussion to G^ v = 
87rGT^ ot we learn that, provided D^T^ = is obeyed, the validity of the tt- and rr- 
component equations (|3.10l a.b) implies the validity of the 99- or y9y9-equation (|3.10cp . In 
the following we replace this latter equation, f|3. l(Jc|) . with D^T^ 1 = as the indepen- 
dent condition. It is satisfied provided the consistency condition (|3.12j) and the ordinary 
continuity equation (|3.5j) are valid. 

Let us summarize. After fixing the form of 9^ and picking a background configu- 
ration G(r), A(r), a system of independent (differential) equations consists of 



These are 5 equations for the 4 unknown functions A(r), B(r), p(r), and p(r). Even though 
the system appears to be overdetermined at first sight it has indeed nontrivial consistent 
solutions. An obvious example is provided by the Classes II and Ilia. From the general 
discussion we know that every solution of the classical Einstein equation induces a solution 
of the modified one. As a consequence, every static isotropic solution of the analogous 



the tt-component (j3.1(Jal) 



the rr-component ()3.10b|) . 



the consistency condition (|3.12j) 



(3.13) 



the continuity equation (|3.5jl . 



the equation of state (|3.4j) . 
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classical equations generates a consistent solution to the above 5 equations. The solution- 
generating algorithm based upon Weyl transformations requires the choice 9^ v = ^ D . 
Based on our experience with cosmological spacetimes [15] we would expect that the 
coupled system has solutions also for other choices of the ^-tensor, 9^ v = for instance. 
(However, an "extreme" choice such as M „ = might be in conflict with the positivity 
properties of T^ v expected for ordinary matter and might necessitate a reinterpretation 
of p and p, see ref. [15].) The weak field approximation of the system (|3.13|) is discussed 
in Appendix 1X1 

Since the Brans-Dicke-type ^-tensor plays a distinguished role [15] we shall adopt 
the following choice for 9^ v in the rest of this paper: 

e ia , = e6™. (3.14) 

Here e is an arbitrary real constant. Particularly interesting values include e = (Class 
I) and e = 1 (Classes II and Ilia). From ()3.14j) with ()2.17|) we obtain a diagonal 9" of 
the type (|3.8|) with the entries 

-^4 (?X> ( 3 - 15a ) 



32vrG A \G , 

p "= e *'=-iki{%) 2 ' < 315b > 

3e 1 /G" x2 



-asciU; ■ (315c) 

We shall see that the essential conclusions of this paper are actually independent of the 
choice for 9^. 

With the choice ()3.14|) we can now write down the differential equations in explicit 
form. We find the tt-component of Einstein's equation 
1 1 A' 



ry 2 ry ^/^2 

= -A-8nGp + 
the rr-component 



2 A' 



rA 2 A 2 



G'\ 1 , N (G'y 1 / G" 
8-3e — + 



G 4:A K \G I A \ G 



(3.16a) 



1 B' 

+ — = -A + 8irGp + 



r 2 A r 2 r AB 



B' 2 

+ 



2AB rA 



G'\ 3e /G" x :! 



G -~A G 1 ■ (al6b) 
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the 99- and ^-component 

A' 



+ 



B' 



A' B' 



B 



/2 



2r A 2 
= -A + 8irGp + 

1 



2rAB AA 2 B 
1 A' 
7A~2A 2 



AAB 2 
B' 



+ 



AA 



and the consistency condition 



A' A' B' 
+ . . n „ + 



B 



,2 



(8-3e) 



B" 



G 



r + 

1 

A 



2AB 
G" 



B" 
2AB 

~G 



(3.16c) 



G 



rA 2 AA 2 B AAB 2 



+ 



3sA' 

Ta 2 



2AB 
1 



G 



+ iXs' 2 - 3E » + ^' 2 - 3£ » 



3e_ f G'G" 
2A V G 2 

2 



+ 



3e 



~G 



2A 

-A' + 8ttG' P 



G 
-- 0. 



(3.17) 



3.3 Selfconsistent cutoff identification 



Up to this point we assumed that G and A are given to us as functions of the 
radial coordinate r. In the RG setting this is not the case, however. The flow equation 
provides us only with a RG trajectory (G(k), A(k), • • • ) into which one must insert a 
cutoff identification k = k(x), or here k — k(r). A detailed discussion of the physical 
cutoff mechanism [15] and its earlier use in the improved solution approach [22] suggests 
an identification k(r) oc l/d(r). Here d(r) is the proper length of a radial curve C, defined 
by dt = d9 = dip = 0, connecting a point with coordinate r to the origin (r = 0). In view 
of the symmetries of the problem the choice of the curve C is certainly the most natural 
one. The more subtle question is: which metric should we use for the computation of the 
proper length d{r)l 

If the renormalization effects are weak and G and A are approximately constant, one 
might pick some background metric 7^, a solution for the case of exact constancy, and 
then compute d(r) = d 7 (r) using this background metric. In the weak field approximation, 
for instance, the Minkowski metric is the natural choice for 7 MJy . In the Classes II and 
Ilia the metric 7^ which generates the actual solution g^ v by a Weyl rescaling can serve 
as a background metric for this purpose. It is clear, however, that the "^-improvement" 
based upon the identification k(r) = ^/d 7 (r) becomes inapplicable when the RG effects are 
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strong. Then g^ u differs significantly from 7^, and since physical distances are measured 
by g^ v one should compute the proper length d(r) = d g (r) from this metric. 

This leads to the complication that in this manner the identification k = k(r) 
becomes functionally dependent on the metric which we do not know beforehand but 
rather would like to compute. In fact, for metrics of the form (j3.1|) the cutoff identification 
of this "^-improvement" involves the function A(r): 



subsection. This results in a system of coupled integro-differential equations for A(r), 
B(r), p(r), and p(r). The input for this system are the fc-dependent "constants" G{k) 
and A(A;). Their r-dependent counterparts are obtained as part of the solution. Given a 
solution for A(r) and B(r), the metric g^ and the corresponding distance function d g (r) 
are known. This in turn fixes the hitherto unknown r-dependence of G and A. In this 
way the cutoff (|3.18j) determines the relation between k and r dynamically. We shall refer 
to it clS cL selfconsistent cutoff identification therefore. 

Solving the system of integro-differential equations 3 appears to be a formidable task 
in general. Therefore it comes as a very pleasant surprise that, even with a selfconsistent 
cutoff, solutions of Classes II and Ilia can be found very efficiently in closed form almost. 
This will be the topic of the next section. 

Before closing this section a comment might be in order about the possibility of 
horizons where the function A(r) turns negative. In ref. [22] a "^-improvement" of the 
Schwarzschild solution has been performed for which A{r) — [l — 2GM/r] 1 is negative 
inside the horizon, and r plays the role of a time coordinate for r < 2GM. In [22] detailed 
physical arguments were given which indicate that, where ds 2 is negative along C, one 
should define the "distance" function as the integral over ds 2 = — ds 2 , rather than ds 2 



3 In principle one can convert the integro-differential equations to differential equations by treating 
d g (r) as an independent function satisfying the additional differential equation d' g (r) = A(r) . 




(3.18) 



Using (|3.18|) we define the r-dependent "constants" G{r) = G(k = A;(r)) and A(r) = 
k(k = k(r)) and insert them into the system of equations (j3.13J) we derived in the previous 
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itself. We adopt the same rule here which amounts to replacing A by \A\ in (|3.18|) : 

r 

d g (r) = J^/\ds 2 \ = Jdr' \l\A{r')\ . (3.19) 
c o 

(See [44] for a related discussion.) However, for the purposes of the present paper the 
question of how to define C inside the Schwarzschild horizon is completely irrelevant. 
We are interested in the physics at r > 2GM where the contribution of the interval 
< r < 2GM to d g (r) is utterly negligible. 

4 Class II vacuum solutions 

In this section we implement the idea of a selfconsistent cutoff identification in 
the framework of the solution-generating Weyl transformations. This restricts us to the 
Classes II and Ilia. Both of them require that T = 0, leading to the "radiation" equation of 
state p = p/3 in the hydrodynamical model. Since this equation of state is not appropriate 
for the physical situations we envisage in the present paper we satisfy T = in the trivial 
way p = p = 0, i. e. we concentrate on vacuum spacetimes. The difference of the Classes II 
and Ilia is that in the first case A(r) vanishes identically, while it is inversely proportional 
to G(r) in the second: A(r) = A G/G(r) ^ 0. Since on the (galactic) scales we are 
interested in the cosmological constant is unlikely to play a decisive role we are not going 
to analyze the case Ilia here. Moreover, one should be aware of the following difficulty. If 
one actually wants to explore the effect of a r-dependent cosmological constant it is not 
enough to have solutions of Class Ilia available. They are indeed realized near a fixed 
point where A(k) oc 1/G(k) oc k 2 so that A(r) oc 1/G(r), but other relevant regimes, a 
"cross over" from the fixed point to ordinary gravity with constant G and A, say, are not 
in Class Ilia and cannot be obtained by the Weyl technique. 

For these reasons we confine the following discussion to vacuum spacetimes of Class 
II, setting A(r) = 0, T^ v = 0, and 6^ u = 6^ D , i.e. e — 1. Then solutions g^ u of the 
modified Einstein equation are obtained by Weyl-transforming solutions 7 M1 , of 

GV(7) = 0. (4.1) 
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According to eq. ()2.21a|) . g^ u = \G(x)/ G\ 7^, or 

ds 2 = [G{x)/G] ds 2 (4.2) 

where ds 2 is the line element of g^ v and ds 2 = ^ ^ dx^ 1 dx u . We implement the symmetry 
requirements by demanding that 7^, too, describes a static isotropic spacetime, and that 
G(x) is annihilated by all 4 Killing vectors, K 11 d^G = 0, implying that it can depend 
on the radial coordinate only. Hence, by Birkhoff's theorem, the only eligible 7^ is the 
Schwarzschild metric 

ds 2 = -f(p) dt 2 + f(p)- 1 dp 2 + p 2 da 2 (4.3a) 

with 

M - 1 " (4.3b) 

P 

For a reason which will become clear in a moment we denoted the radial coordinate p 
rather than r. (There should be no confusion with an energy density.) Note that the factor 
G in (|4~3bjl arises because the coefficient of the 1/p-term is a free constant of integration, 
usually written 2GM, and to be interpreted as 2GM in the present context; G does not 
appear in the field equation 6^(7) = 0. Minkowski space, the special case M = 0, can 
be used as well for generating solutions. 

If G was known as a function of p, the line element we are after could be obtained 
as ds 2 = [G(p)/ G] ds 2 simply. However, as we discussed in Subsection 13.31 what the 
RG analysis provides us with is the function G(k) only. In the following we employ the 
self consistent, hence g^-dependent, cutoff identification 

Hp) = i/d g { P ) (4.4) 

in order to convert G(k) to G(p) = G[k = k(p)). 

Given the function G(k), it will prove convenient to define an associated function 
W(y) by 

W{y) = G{k = i/y)/G (4.5) 

with £ an arbitrary fixed constant of order unity whose precise value will not matter. 
Obviously W (y) with y = d g (p) is the Weyl factor relating the two metrics: 

ds 2 = W(d g {p)) [-/(p)dt 2 + /(p)- 1 dp 2 + p 2 d ( 7 2 l. (4.6) 
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Eq. ()4.6|) is an implicit equation for g^: on the RHS of ()4.6|) it makes its appearance in 
d g (p), on the LHS in ds 2 . Solving this equation is tantamount to finding d g as a function 
of p. In fact, once we know the function p i— > d g (p) we can combine it with y ^ W(y), 
known from the renormalization group, and define 

w{p)=W(d g {p)). (4.7) 

By definition, to is a function of p, while the natural argument of W is y = d g (p). We 
have accomplished our task when we know w , for then 

ds 2 = w{p) [-f(p) dt 2 + f{p)- 1 dp 2 + p 2 da 2 ] (4.8) 

is an explicit formula for the line element of g^ u . 

How can we get a handle on d g (p)7 Eq. (|4.6|) tells us that, if dt = d6 = dip = 0, 

\ds 2 \=W(d g (p))\f(p)\~ 1 dp 2 (4.9) 

since G and W are positive throughout. The total proper length of C obtains by integrating 

p 

d g (p) = [dp' ^W(d g {p'))\f{p>)\- 1 . (4.10) 



o 



Since W and / are known functions (J4.1U)) is an integral equation for d g (p). Differentiating 
it yields the equivalent differential equation 



d 



d a (p) = yjw(d g (p))\f(p)\- 1 , d g (0) = 0. (4.11) 



dp 

All we have to do in order to implement the selfconsistent cutoff is to solve this differential 
equation for d g (p). Then, by (|4.7j) . we know w(p) and, by (|4.8|) . the solution g^ u . 

Eq. (|4.11j) is easily solved by a separation of variables. After computing the integrals 

y 

I w (y) = [ (4.12) 
W(y') 



o 

p 



lj(p) = / -M= (4.13) 

7(p0| 
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it remains to solve the equation Iw(d g (p)) = If(p) for d g (p): 

d g (p) = I w 1 (l f (p)). (4.14) 

Note that Iw depends only on the RG trajectory, while If requires information 
about 7 Mi , only. With (|4.3b|) the latter integral is elementary. For p > 2GM: 



I f {p) = %GM + 2GM In 



pj (2GM) + Jp/ (2GM) - 1 



(4.15) 



+ y/p(p- 2GM) . 
The asymptotic behavior for p large is 

I f ( p ) = p + G M In (p/p ) + const + ( (G M) 2 /p\ (4.16) 

with an irrelevant constant p Q . 

Our final result, the line element f)4.8|) . while describing a static isotropic spacetime 
is not in the standard "A£>-form" ()3.1|) . For some applications it is advantageous to 
analyze it in the (t, p, 8, (/?)-system of coordinates as it stands. For others it is better to 
convert it to the AB-form by an appropriate change of coordinates. It is sufficient to 
replace p with the new radial variable r defined as 

r(p)=p^(p). (4.17) 

Assuming that the relationship r = r(p) can be inverted to obtain p = p(r), it is straight- 
forward to rewrite ()4.8j) in the (t, r, 6, (^)-system: 

ds 2 = -wfdt 2 + +r 2 da 2 . (4.18) 

/ [l + ld p \nw] 2 

The components of the metric ()4. 18|) are to be evaluated at p = p(r). 

Summarizing the results of this section we can say that vacuum solutions of Class 
II, with a selfconsistent cutoff identification, can be obtained by means of the following 
algorithm. 

Step 1: Select a trajectory G{k) and write down the associated function W{y) of ()4.5|) . 
Step 2: Perform the integral ()4.12j) to obtain Iy^ and find the inverse function J^) 1 . 
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Step 3: Write down d g (p) according to (|4.14p . Use it to construct w(p) = W(d g (p)). 



After having completed Step 3 we know the final result for the metric in the 
(t, p, 9, (^)-system, eq. (|4.8|) . Only at this point we learn what is G as a function of p. 
Since G(p) = G(k = £/d g {p)) = GW(d g (p)) we have 



G{p)=Gw{p). 



(4.19) 



This is the position dependence of Newton's constant which, by means of the selfconsistent 
cutoff identification, gets associated to the trajectory G = G(k) in a dynamical way. As 
we shall see, the interrelation between G{k) and G(p) can be highly nontrivial and could 
not be guessed beforehand. 

If one needs the metric in standard form one continues with 



Step 4: Define r(p) = p \Jw(p) and find the inverse function p = p(i 
Step 5: Read off the components of the v4.B-metric from ([4. 18)1 . 



A(r) = f( P y 1 [l + fohxwQ*) 
B(r) = w(p(r)) f(p(r)), 



and express Newton's constant function of r: 



p=p(r) 



(4.20a) 
(4.20b) 



G(r) = Gw(p(r)) =G 



p{r) 



(4.21) 



It can be verified explicitly that, for every given RG trajectory G = G(k), eqs. 
f)4.20l a.b) along with ()4.21|) do indeed constitute a solution to all components of Einstein's 
equation, eqs. (|3.16l a.b.c) and the consistency condition (|3.17|) for the case T^ v = 0, A = 0, 
and O^u = . The calculation is quite lengthy and tedious. It demonstrates explicitly 
that, as it was to be expected [15], the algorithm works only for e — 1, i.e. for the 
Brans-Dicke ^-tensor. 
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5 The IR fixed point scenario 

5.1 Scale dependence of G 

In this section we apply the algorithm for obtaining Class II vacuum solutions de- 
veloped in the previous section to the scale dependence 

G(k) = G + g l */k 2 . (5.1) 

As we discussed already in the Introduction, it interpolates between classical gravity 
G(k) = G for large k, and the fixed point behavior G(k) oc 1/k 2 in the IR (k — ► 0), the 
transition taking place near the scale k 2 r = / G at which the two terms on the RHS of 
()5.1|) are exactly equal. The choice (j5.1|) is motivated by its phenomenological success in 
cosmology [17, 18]. 



5.2 Class II vacuum solutions with selfconsistent cutoff identi- 
fication 

Let us now follow the steps listed at the end of Section HJ According to ()4.5j) . the 
trajectory (|5.1|) amounts to 

L 2 

with 



W(y) = l + K 2 y 2 = l + ^- (5.2) 



1 o IR 

* 2 - £ 'fir- < 6 - 3 > 

Here L is the proper distance d g at which the fixed point behavior sets in, and k = L^ 1 
denotes the corresponding mass scale. For the function ()5.2j) the integral Iw of ()4.12|) can 
be performed in closed form: 



I\v(y) = K 1 arsinh (ny) = k 1 In 



(5.4) 



The next step consists in finding d g (p) from the implicit equation I w (d g (p)) = If{p)- 
For (|5.4JI this is easily done: 

d g (p) = tC l sinh(/« If(p)) ■ (5-5) 
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Recalling eq. (|4.15|) for If, we see that d g (p) is known in completely explicit form now. 
By fl4~7|) . it implies the Weyl factor w(p) = W{d g {p)) = 1 + [nd g (p)\ 2 1 or, 

w(p) = cos1i 2 (k If(p))- (5.6) 

Thus we are in a position to write down the final result for the metric and the 
position dependence of G. From (|4.8J) and ()4.19|) we obtain, respectively, 

ds 2 = cosh 2 (K I f (p)) \-f(p) dt 2 + f(p)~ l dp 2 + p 2 da 2 ] (5.7a) 

and 

G{p) = G cosh 2 (k If{p)) . (5.7b) 

Actually (j5.7|) describes a one-parameter family of solutions labeled by the parameter M 
which enters via f(p) = 1 — 2GM/p. 

Our construction includes M = as a special case. Here If(p) = p entails the 
simplifications 

d g (p) = kT x sinh (/t p) , (5.8a) 

w(p) = cosh 2 (k p) , (5.8b) 

G(p) = G cosh 2 (Kp), (5.8c) 

ds 2 = cosh 2 (k p) [-dt 2 + dp 2 + p 2 da 2 ] . (5.8d) 

Eqs. ()5.8|) describe a curved, but conformally flat, stationary and isotropic spacetime. 
At the origin (p = 0), Newton's constant has the value G, and at distances p > L it 
starts growing very rapidly in an exponential fashion. Since the conformal factor w grows 
likewise without any bound, the spacetime is not asymptotically flat. 

Eq. ()5.8a|) shows that for p > L the proper length d g differs significantly from the 
coordinate "distance" p. This explains the a priori completely unexpected p-dependence 
of G, eq. (J5.8c|) . It is the result of the selfconsistent cutoff identification which determines 
k = k(p) dynamically. With a naive cutoff employing the distance function <i 7 (p) = p of 
the Minkowski metric 7 Miy we would obtain G(p) = G [1 + k 2 p 2 ] instead. This coincides 
with (|5.8cj) for p <^ L, but differs substantially for p^> L. In the latter regime it would 
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make no sense to use the naive cutoff because the renormalization effects are strong, w(p) 
is large, and rf 7 (p) no longer has anything to do with the physical distance. 

For M > the solutions ()5.7|) have a similar structure, with an additional distortion 
due to the curvature of the Schwarzschild spacetime. It becomes strong for p — > 0, but in 
this regime our model is anyhow inapplicable because the trajectory (|5.1|) does not "know" 
about the UV running of G(k). In the following we shall analyze the solution ()5.7|) only far 
away from the horizon (p ^> 2GM) . Typical question about black hole physics (structure 
of horizons, nature of the central singularity) must involve a more complete description 
of G(k) and will be addressed elsewhere. 

Let us return to the general solution for M > 0. In (|5.7|) it is written in the 
(t, p, 9, (^)-system of coordinates. Introducing the new radial coordinate r = p yVw(p) we 
can transform it to the standard AB-form, but unfortunately it is not possible to obtain 
a closed formula for p as a function of r. Eqs. (I4.19J1 . ()4.20|) yield 

-2 



A(r) 



B(r) 




G(r) = G 



— I T 



p{r) 



(5.9a) 

(5.9b) 
(5.9c) 



The function p(r) is defined by the transcendental equation 

r = p(r) cosh nlf(p(r)) . (5.9d) 

In order to get a first understanding of the physics behind ()5.9|) we are now going 
to analyze the regimes 2GM <C p <C L and 2GM <C L p in turn. 



5.3 The modified Newtonian regime: 2 G M p <^ L 

We consider the portion of spacetime where p is much larger than its value at the 
Schwarzschild horizon, 2GM, and at the same time much smaller than the character- 
istic scale L. In the application to galaxies which we have in mind, 2GM is basically 
negligible compared to L so that p can vary over many orders of magnitude without vi- 
olating 2GM <C p L. Because here 2GM/p <C 1, as in the Newtonian limit of the 



28 



Schwarzschild metric, we call this regime "modified Newtonian" . What we are interested 
in are the modifications due to the running of G. 

Let us expand the solution ()5.9|) . with If(p) given by ()4.16|) . for small values of 



G M GM 



Inp, (9 2 «1. 



(5.10) 

P p ■ \LJ 

We retain only terms linear in the quantities (|5.1()jl and neglect higher powers and cross 
terms. In leading order one has p{r) = r [l — | k 2 r 2 + ■ ■ ■ ] , and ()5.9|) yields 



A{r) = 1 + 

B(r) = 1 - 
G{r) = G 



2GM 
r 

2GM 

r 



2n 2 r 2 , 



, 2 2 
+ K T , 



1 + k 2 r 2 



(5.11a) 

(5.11b) 
(5.11c) 



Since p, d g , and r differ only slightly, r is basically the same as the physical distance so 
that (I5~TTJ) is valid for 2GM <r<L. 

In order to interpret the above result and similar metrics which we shall encounter 
later on let us recall a well-known fact about the Newtonian limit [43]. Assume we are 
given a metric of the form 



ds 2 



1 + 2 0) dt 2 + (5 i:i + hij) dx* dx j , 



(5.12) 



valid to first order in and hij. The timelike geodesies of test particles are obtained from 
the mechanical action S m = —m fy/—ds 2 = —m Jdt [1 + 20— (5ij + hij) x l x^] 1 ^ 2 where 
x l = dx l /dt. In the Newton limit one assumes that the particle moves at nonrelativistic 
velocities v 2 = 5{j x % i J < 1, and that 0, hij, and v 2 are small of the same order. Expand- 
ing S m to first order in those quantities leads to S m = Jdt {|mt> 2 — m0}. This shows 
that the function in 1)5.12)1 is indeed the Newtonian potential, and that the Newton 
limit is insensitive to the spatial components hij since hij x % x- 7 is a second-order term. 

By (j5.11j) . the Newtonian limit of our static isotropic spacetimes is described by the 
line element 



ds 2 



2GM 



i 2 2 

+ k r 



dt 2 + 



1 + 



2GM 



2 k 2 r 2 



dt 2 + r 2 da 2 . 



(5.13) 
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According to the above discussion, it implies the following modified Newtonian ("mN") 
potential: 

mN (r) = -^ + i*; 2 r 2 (5.14) 

While m N shows the familiar 1/r-behavior at small distances, it becomes positive and 
grows proportional to r 2 at large distances. The term oc r 2 is a direct consequence, and 
in fact the only one in the Newtonian limit, of the IR fixed point running of G. It has 
a rather dramatic consequence: at sufficiently large distances, the force two masses exert 
on each other grows linearly with their distance. This behavior is reminiscent of the 
confinement forces among quarks in QCD. 

What is the range of validity of the (modified) Newtonian approximation? Is it really 
applicable to a regime of radii where the r 2 -behavior is already fully developed, overriding 
the classical 1/r-term in (j5.14j) ? The answer is indeed in the affirmative. The reason is 
the simple but very important observation that the nonclassical ("nc") behavior already 
sets in at a radius r nc which is much smaller than the scale L (where the approximation 
breaks down). 

Let us define r nc by the condition that the two terms on the RHS of (j5.14j) are 
exactly equal there: GM/r nc = /t 2 r 2 c /2. For r < r nc the classical term is dominant, 
while the r 2 -term prevails for r > r nc . Writing r s = 2GM for the Schwarzschild radius, 
the transition to the nonclassical r 2 -regime takes place at 

r nc = (r s L 2 f 3 . (5.15) 

As a consequence, the ratio r nc /L is 

In fact, the modified Newtonian limit applies when rs C r < L, so that rs L 
in particular. Thus we conclude that there is a wide range of distances, r nc < r <C 
L, where the nonclassical behavior dominates m N without being in conflict with the 
approximations made. 

The natural emergence of a "hierarchy" r nc C L is a very important feature. For 
instance, when we apply this theory to galaxies later on r nc should be a typical galactic 
scale while L is of the order of the Hubble radius. 
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5.4 The fixed point regime: 2GM«I« J o 



In this section we analyze the solution (|5.7j) or, equivalently, ()5.9|) at asymptotically 
large distances from the center. As the Schwarzschild metric is asymptotically flat, we 
can approximate I/(p) ~ p for p — > oo. As a result, all terms containing M drop out, 
and it is sufficient to study the k p — > oo-behavior of the M = 0-solution, eqs. ()5.8|) . In 
this limit we can easily transform the metric to the standard AB-foim. Since cosh (k p) « 
exp (ftp) /2 for «p — > oo, the implicit equation for p = p(r), eq. (|5.9dp . becomes 2r = 
p(r) exp[/cp(r)]. It can be solved in terms of the Lambert W-function 4 Wl: 



p(r) = k 1 W L (2/cr) , Kr»l. 
Hence the asymptotic form of the solution, in the (t, r, 6 1 , <p)-system, reads 

A(r) = [l + W L (2nr)]~ 2 , 

B(r) 



k r 



1 2 



G{r) = G 



K T 



(5.17) 

(5.18a) 
(5.18b) 

(5.18c) 



For x — > oo there exists an asymptotic expansion of Wl in terms of iterated logs [46]: 



,aw n i in In In a; 

VVl(^) = In a: — minx H — h 

mx 



(5.19) 



Therefore, at very large distances, B{r) grows proportional to r 2 / In 2 (2 nr). The asymp- 
totic spacetime geometry is strongly curved therefore. This fairly exotic regime will not 
be relevant in the context of galaxy haloes. 



5.5 Geodesies 

In order to gain a first understanding of the dynamics implied by the metric of 
eqs. (|5.9|) it is instructive to investigate its timelike geodesies, the trajectories of massive 
test particles. Here we concentrate on periodic circular orbits at constant radius r or p. 

4 By definition [46], the Lambert W-function satisfies Wl(x) exp [W L (:r)] = x. In our notation, Wl 
stands for its real branch analytic at x = 0, denoted W or Wo in ref. [46]. 
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In the Newtonian limit the analysis is very simple. The modified Newtonian poten- 
tial m N allows for circular orbits at any radius r. On those orbits the centrifugal force 
exactly balances the gravitational attraction provided the velocity has the time indepen- 
dent value 

v \r) = r <f)' mN (r) = — + k 2 r 2 . (5.20) 

At small radii r < r nc , the classical term G M/r dominates, yielding v(r) oc 1/ ' y/r for the 
velocity and Kepler's T(r) oc r 3//2 for the period of the orbit, T(r) = 2nr/v(r). Instead, 
for r > r nc , the nonclassical term becomes dominant, giving rise to a linearly increasing 
velocity v(r) — kt and a period T(r) = 2tt/k = const which is independent of r. 

Note that there exist periodic orbits even if M = 0. In this case they owe their 
existence entirely to the r-dependence of G since the solution-generating metric 7^ with 
M = (Minkowski metric) does not give rise to periodic orbits. 

Going beyond the Newtonian limit, the analysis is most easily performed using the 
radial coordinate p rather than r. In Appendix El we derive the geodesic equation for a 
metric of the type (J4.8)) with arbitrary functions w (p) and f(p). We find that it admits a 
circular orbit at coordinate radius p provided the quantity 



1 1 - , (5.21) 

\2n P J f' + ^f 

is positive. If so, T(p) is the period of the orbit, measured in units of the coordinate time 
t. 

In the IR fixed point scenario the Weyl factor w of the exact metric is given by (|5.6J) 
and / is of the Schwarzschild form (|4.3bjl . Inserting these functions into the equation 
(j5.21j) one finds that its RHS is indeed positive for any p > 0. Hence the exact solution 
admits a circular orbit at every radius. 

A simple, closed-form expression for the period T(p) can be obtained for M = 
only. In this case, according to (|5.8|) . / = 1 and w(p) = cosh 2 (up) which yields when 
inserted into (j5.21|) : 



T(p) = 2np 



coth (k p) 



1/2 



(5.22) 



(np) 

The M = 0-spacetime ()5.8d|) is regular at p — 0, and the Weyl factor equals unity 
there. Therefore the coordinate time t equals the time indicated by a clock which is at 
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rest at p = 0. Correspondingly, T(p) is the period as measured by this specific clock. We 
saw already that for p — > oo the spacetime is strongly curved and that the proper distance 
to the origin, d g (p) = y, is exponentially large compared to the coordinate distance p. In 
fact, by (|5.8aj) . ny — sinh (k p). It is more physical, therefore, to express T in terms of y 
rather than p: 



2n 



T(y) = - 



K 



arsinh^ (k y) + arsinh (« y) \l 1 + (ny) 2 



1/2 



(5.23) 



At small distances y <C = L, this function has a power series expansion T(y) = 
2n k~ 1 [l + | (ny) 2 + 0((k?/) 4 )] which, for y — > 0, reproduces the constant period found 
in the Newtonian limit: T »s 27r/t~ 1 . Conversely, at extremely large distances y ^> «~ , 
one finds a growth of T(y) which is approximately logarithmic: 

T(y)& — ln(2«y). (5.24) 

As we mentioned in Subsection 15.41 already, all terms containing the parameter M 
are negligible in the fixed point regime 2GM <C L p. As a consequence, for y — > oo, 
also the solutions with M^O give rise to the asymptotic behavior (|5.24jl . 



5.6 Cosmological fixed point running at galactic scales? 

In this subsection we address the question of whether the non-Gaussian IR fixed 
point which was invoked in refs. [17, 18] in order to describe the cosmology of the late 
Universe might also be responsible for the non-Keplerian behavior of the observed galaxy 
rotation curves v(r). The model trajectory (|5.1|) describes "ordinary" gravity with G = 
G = const at small distances, i. e. at the scale of laboratory or solar system measurement, 
and the IR fixed point running at cosmological scales. The typical scale of galaxies, 1 kpc 
say, is in between those two regimes. The intriguing question is whether the unexpected 
features of their rotation curves can be interpreted as due to the onset of the cosmological 
fixed point running. If so, the parameter g^ 1 appearing in our present ansatz ()5.1|) must 
have the same numerical value as in cosmology since one and the same RG fixed point 
would govern both cosmology and galaxies. 

In the extended IR fixed point model proposed in ref. [18] the fixed point era proper, 
with G(k) oc l/k 2 , is preceded by a standard Friedmann- Robertson- Walker epoch with 
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G = const, the transition taking place at a redshift z tI . In this model, can be expressed 
in terms of z tr and the present Hubble radius /h(^o) 5; 




(5.25) 



Here Zpi = V G is the Planck length defined in terms of the Newton constant G measured in 
the laboratory. The constant £ relates k to the cosmological cutoff 1/t and may be assumed 
of order unity therefore. Moreover, confronting the model with various cosmological 
standard tests [18] it turns out that z tr cannot be much larger than unity. This information 
allows us to estimate the scale L = 1/k. Using (|5.25|) in (|5.3j) we obtain 

L = ^i(l + z tr )- 3/4 Z H (to). (5.26) 

As an order of magnitude estimate we can say that L equals approximately the present 
Hubble radius since the prefactor of /h(^o) on the RHS of ()5.2(ij) . while not known exactly, 
is certainly of order unity. For simplicity we shall assume in the following that L = Zh(^o); 
hence n = H is presicely the present Hubble parameter, H = 1//h(^o)- 

This result has an important implication when it comes to applying our solution of 
the modified Einstein equation to the galaxy problem. We had analyzed it in the modified 
Newtonian regime p <C L and for p 3> L where the fixed point behavior is already fully 
developed. With L equal to the Hubble radius the latter regime is of no relevance for the 
structure of galaxies. We conclude that the portion of the spacetime which is interesting 
regarding the problem of the rotation curves lies entirely within the Newtonian regime 
p < L. 

In this regime the modified potential of a point mass M is given by (j5.14|) ; its 
nonclassical term grows proportional to r 2 , implying a linearly increasing force and 

v(r) = ztr = H r. (5-27) 

We mentioned already that the scale L may not be confused with the distance r nc 
at which the nonclassical force oc r starts dominating Newton's l/r 2 -force. Eq. (J5.16J) 
tells us that r nc / L = (2GMHoj . To be explicit, let us consider a point mass with 

5 Eq. lj!T25j> obtains by inserting cqs. (17) and (51) of ref. [18] into eq. (29) of [18]. 
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M = 1O 11 M , the mass of a typical galaxy. Using Hq = 70 km sec 1 Mpc 1 we find the 
ratio 

r nc /L « 1.3 • KT 4 (5.28) 

which amounts to 

r nc ^ 560kpc. (5.29) 

Several comments are in order here. First of all it is gratifying to see that our 
model is indeed capable of producing a hierarchy between the cosmological scale L and 
the (would-be) galactic scale r nc . However, quantitatively this hierarchy is still too small 
by about 2 orders of magnitude to be phenomenologically acceptable; rather than ()5.29|) . 
we would like to find a r nc of the order 1 kpc, which is roughly the distance at which the 
plateau behavior in the rotation curve of a typical spiral galaxy sets in. In a real galaxy 
the luminous matter has a complex spatial distribution while above we were studying the 
field of a point mass. However, it appears unlikely that this is the origin of the missing 2 
orders of magnitude. We rather take them as a hint indicating that, while qualitatively 
pointing in the right direction, the RG trajectory (|5.1|) is probably not very realistic at 
the kiloparsec scale. 

An even stronger argument for rejecting this trajectory comes from the rotation 
curve it gives rise to, eq. (|5.27|) . A linearly rising v(r) for a point mass can hardly be 
reconciled with the almost constant velocities observed in galaxy haloes. Moreover, the 
predicted velocities are substantially too small. In fact, (J5.27j) happens to have the form of 
Hubble's law, and Hq r is several orders of magnitude smaller than the measured plateau 
velocities of a few hundred km/sec, of course. 

The results of this section suggest that, at galactic scales, G(k) should have a milder 
growth for decreasing k than it has in the IR fixed point model. A realistic RG trajectory 
should entail a m N which increases with r less steeply and makes v (r) approximately 
constant at large distances. At the same time it must produce higher orbital velocities. 
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6 Power law running 



6.1 The RG trajectory 



The discussion at the end of the previous section has indicated that an infrared 
growth of G(k) proportional to 1/k 2 is probably too strong at the scale of galaxies. This 
motivates us to investigate the consequences of the following model-trajectory with a 
milder fc-dependence of G(k). We assume that G(k) is approximately constant for large 
k, and that it runs according to the power law G(k) oc l/k q at small k. Here q is 
an arbitrary positive constant. Applying the model to galaxies we shall find that the 
relevant range of g-values is < q <C 1 so that the variation of G(k) is indeed much 
weaker than in the IR fixed point scenario. 

For the time being we allow for arbitrary q > and make the following ansatz: 



Here G and the transition scale k tr are free constants. Note that G(k) is continuous 
but not differentiable at k tr . Recalling that the dimensionful Newton constant evolves 
according to k dj-G = 7/n G, we see that the model trajectory (j6.1|) amounts to a constant 
anomalous dimension rj^ = —q for k < k tr and to rj^ = at scales k > k tr . A more 
realistic RG trajectory would interpolate smoothly between 7/n = — q and 7/n = 0, but for 
our present discussion the idealization of a jump at some k tT is good enough. Our main 
conclusions will be insensitive to the intermediate regime interpolating between standard 
gravity and the power law regime. 

6.2 Class II vacuum solutions 

Next we apply the Weyl-algorithm of Section 0] to the G(k) of eq. (|6.1|) in order to 
find the vacuum solutions of Class II with a selfconsistent cutoff identification. The first 
step consists in identifying the function W(y). From (J6.1j) we obtain 





for k > k s 



tr 



tr 



(6.1) 





for y < L 



(6.2) 
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For ktr and £ fixed, we introduce the length scale L and the mass scale k = L 1 by 



k tr = — = £ k. 



(6.3) 



(The quantities k and L are analogous to those of the previous section but not mathemat- 
ically identical.) The cutoff identification k(p) = i/d g {p) provides us with the p- value p tr 
which corresponds to the transition scale k tT ; from k(p tT ) = k tr we obtain the condition 
d g {ptr) = L. Since £ is of order unity, n is essentially the same as k tr and L is its inverse. 
When measured with the physical metric g^ ui the length scale L is the proper length 
corresponding to the coordinate value p tT . 

For W(y) from ()6.2|) the integral Iw of (j4.12|) can be done easily: 



Iw{y) 




(1-9/2) 



(!) 



l-q/2 



for y < L 
for y > L 



(6.4) 



According to ()4.14|) we obtain d g (p) = /^(//(p)) by inverting Iw- The result reads 



d g {p) 



(1-1) 



i tip) 



+ - 



2-g 



for p < p tr 
for p > p tr 



(6.5) 



Since d g = If for p < p tr , the condition dg(p tr ) = L defining p tr simplifies: 



J /(Pt r/ 



L. 



(6.6) 



Finally we can obtain the Weyl factor w(p) = W[d g (p)) by inserting (|6.5|) into (|6.2j) : 



w(p) 



(1-f) «//(p) + 



2-9 



for p < ptr 
for p > ptr 



(6.7) 



With the derivation of w(p) we have completed the construction of a one-parameter 
family of solutions to the modified Einstein equation with = and A = 0. The 
solutions are labeled by the mass M. In the (t, p, 8, (^)-system of coordinates, the metric 
is given by ()4.8|) with ()6.7|) and the position dependence of Newton's constant by G(p) = 
Gw{p). 
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6.3 Rotation curves 



In this subsection we study timelike geodesies of the metric with the Weyl factor 
w(p) of ([6.7)1 . Again we consider circular orbits only, being particularly interested in their 
rotation curves. As for the very existence of circular orbits, we showed in Appendix |E] 
that for this class of metrics there exists an orbit at constant p provided the RHS of 

2 



(T(p)Y_ 1+Y(p) 
\2«p) (GM/p)+f( P y 



(GM/p)+f(p)Y( P y 

with Y(p) = \ pw'(p)/w(p), is positive at this value of p. For p < p tT we have Y = and 
the situation is the same as with the Schwarzschild metric: If M > there are circular 
orbits, satisfying Kepler's law (T/2n) 2 = p 3 / (GM), and there are none if M = 0. In the 
more interesting region p > p tr the quantity Y(p) is given by 

Y{p) = t qKp , ; ■ (6.9) 

[(2-q)Kl,{p)+q] y/M 

It can be shown that the RHS of ()6.8|) with (|6.9J) is a finite positive quantity, both if 
M > and M = 0. Hence there exist circular orbits at all radii p > p tr and their periods 
are given by T(p). Such orbits exist even if, for M = 0, the construction of the solution 
g^u is based upon the Minkowski metric which has no circular orbits. In this case the 
gravitational attraction is entirely due to the gradient of Newton's constant. 

It is instructive to discuss the following special cases which can be dealt with ana- 
lytically. 

(i) Asymptotically large radii 

Let us focus on orbits with very large radii p — > oo. In this limit, f(p) — > 1 and 
Ij(p)/p — > 1 which implies Y(p) — > qj (2 — q). Using this asymptotic behavior in ()6.8|) 
we see that 



lim 

p— >oo 



T{ P y 2 



2 

6.10 

q 



2vrp 

This is a very remarkable result. In fact, let us tentatively assume that for p — > 
oo the curvature of spacetime is small and that the test particles on their orbits move 
with velocities much smaller than the speed of light. In this case a quasi-Newtonian 
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interpretation is possible which neglects the difference between proper and coordinate 
time and likewise between the proper circumference of the orbit and the euclidean 2-7rp. 
Therefore (27rp) /T(p) = v(p) is "the" velocity of the test particle. In this quasi-Newtonian 
limit, eq. (|6.10j) tells us that for p — ► oo the velocity v(p) approaches a finite, p- independent 
value Voq = limv(p): 

v 2 oa = \q. (6.11) 

Thus we may conclude that, provided the quasi- Newtonian interpretation is applica- 
ble, the scale dependence G{k) oc l/k q does indeed give rise to a rotation curve v(p) with 
a plateau behavior asymptotically. The plateau value = \fqj2 is a direct measure of 
the exponent appearing in the power law. It does not depend on M. 

Clearly it is very tempting at this point to speculate that the observed plateau 
behavior in typical rotation curves of spiral galaxies is due to a power law running of G. 
For the Newtonian picture to apply it is necessary that which implies that the 

exponent q must be much smaller than unity: g < 1. It is gratifying to see that this 
condition is indeed met by the "test particles" in the galactic haloes whose velocity was 
measured [30]. Typical values of range from 100 to 300 km/sec so that, in units of the 
speed of light, Voo ~ 10~ 3 . Hence, ignoring factors of order unity, we find that the data 
require an exponent of the order 

q « 10~ 6 . (6.12) 

Later on in Subsection 16.41 we shall see that also the second prerequisite of the Newtonian 
interpretation, small spacetime curvature, is satisfied for q as small as 10~ 6 . 



(ii) The case M = 

For non- asymptotic radii p the rotation curve v(p) can be obtained in closed form 
only if M = 0. In this special case we have f(p) = 1, If(p) = p, and ()6.8|) becomes 




( P > Ptr) (6.13) 



The RHS of this expression is manifestly positive, showing that there exist circular orbits 
at any p > p tI , with period T(p) = 2irp a/ (2/q) + 1/ (k~p)- Interpreting this result a la 
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Newton yields the velocity 



For p — ► oo it approaches the plateau value v OQ = q/2 from below. 

(iii) First order approximation in q and GM/p 

In order to obtain an analytic, albeit approximate expression for v(p) in the generic 
case M ^ we shall assume that q <C 1 and GM/p <C 1, the second condition meaning 
that p is much larger than the Schwarzschild radius. In the halo of a typical galaxy both 
conditions are well satisfied. In evaluating ()6.8j) we retain only terms of first order in q 
and GM / p and neglect all higher powers and cross terms. The result is then 

v 2 (p) = —+ q - (p>Pt r ) (6.15) 

p A 

This v 2 (p) approaches the plateau value q/2 from above. 

If p is small, the first term on the RHS of ()6.15|) dominates, while the "nonclassical" 
contribution q/2 prevails for large values of p. We introduce a radius p nc which marks the 
onset of the nonclassical behavior. We define p nc to be the radius where the two terms on 
the RHS of ()6.15|) have equal magnitude: GM/p nc = q/2. Obviously the "new physics" 
sets in at 

Pnc=ps/q (6.16) 

where ps = 2GM is the Schwarzschild radius pertaining to 7^. 

Since (16.15)) holds for p > p tT only, the estimate ()6.16|) is valid only provided p nc > 
p tT - However, it is remarkable that, as long as p tr < p nc , the value of p nc is completely 
independent of p tr and the scale L = d g (p tr ). Thus we can say that, at least when 
the quasi-Newtonian interpretation makes sense, the distance scale where the galactic 
plateau behavior sets in, typically 1 kpc in Nature, is quite unrelated to the "cross-over 
scale" k = L^ 1 where the power law running of G begins. 

6.4 Properties of the metric 

In order to fully establish the picture sketched in the previous subsection we must 
justify the validity of the Newtonian approximation in the "halo". This is most easily 
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done using the standard ("AB") form of the metric. At the end of Subsection 16.21 we 
found the Class II solutions pertaining to G(k) oc l/k q in the p-coordinate system, the 
Weyl factor which determines both and G(p) being given by eq. ([6.7)1 . In general it is 
not possible to express this solution in the r-coordinate system in closed form. However, 
there are two instructive limiting cases where this is indeed possible: asymptotically large 
radii and the first order approximation in q and GM/r. Let us discuss those cases in turn. 

(i) Asymptotically large radii 

Retaining for p — > oo only the leading-order terms one has f(p) = 1 + • • • , I/(p)/p — 
1 + • • • , and from (|6.7|) it follows that 

w(p)=[(l-q/2)Kp] 29/{2 - q) . (6.17) 

Using this approximation for the Weyl factor, r = p \Jw(p) yields the following relation- 
ship between the radial coordinates r and p: 

p= [(l-g/2) K]- q/2 r^ 2 . (6.18) 

Eqs. fdH) and (j6T8j) imply w(p(r)) = [(1 - q/2) Kr] q and pd p \nw = 2 q/ (2 - q). As 
a result, ()4.18|) gives rise to the following metric in the (t, r, 9, <p)-system of coordinates, 
valid for r — > oo: 

ds 2 = -[(1 - q/2)nr] q dt 2 + (1 - q/2) 2 dr 2 + r 2 da 2 . (6.19a) 

This solution describes the spacetime geometry due to the r-dependent Newton constant 

G(r) = G [{1 -q/2)K] q r q . (6.19b) 

Note that in deriving the asymptotic forms ()6.19l a.b) we did not assume that q is 
small, they are valid for all q G (0, 2). (As it should be, r increases monotonically with p 
in this interval.) 

The metric (I6.19a|) is the universal r — > oo-asymptotics of all Class II vacuum space- 
times, i. e. it does not depend on the parameter M. This geometry is not asymptotically 
flat which, in fact, is not too surprising since G(r) cx r q grows unboundedly for r — » oo. 
The spatial part of (|6.19aj) is the metric of a cone; the singular tip of the cone (r = 0) 
lies outside the domain of validity of this metric, however. 
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It is interesting to consider the asymptotic solution (|6.19|) for small values of q: 
< (7 <C 1. Expanding in powers of q and neglecting terms of order q 2 and higher yields 



ds 2 = - [l + q ui(k r)] dt 2 + (1 - q) dr 2 + r 2 da 2 , (6.20a) 
G(r) = G[l + q\n(Kr)]. (6.20b) 

One must be careful about the range of validity of (|6.20p . In order to obtain this ap- 
proximation we used the expansion (nr) 9 = exp[g ln(«r)] = 1 + q m(« r) + 0(q 2 ln 2 (ftr)) 
so that the actual expansion parameter is q ln(/tr) and not q itself. This approximation 
is reliable provided q ln(nr) <C 1 or, equivalently, r -C L exp(l/g). Since on the other 
hand p, or what is the same if q is small, r, should be much larger than L, we see that 
()6.20|) is valid in the window 

L < r < L exp(l/g). (6.21) 

For galaxies the upper bound on r is of the order exp(10 6 ) L. Unless L is anomalously 
small it is probably well satisfied for all practical purposes so that we may regard ()6.20j) 
as the asymptotic form of the solution. 



(ii) First order approximation in q and GM/r 

Next we investigate the solution (|6.7|) for all p > p tT , including non-asymptotic radii. 
We expand in q 1 and G Mj p 1 up to the first order discarding higher powers and 
products of these quantities, as well as terms of order q 2 In 2 (up). In this approximation, 



w(p) — 1 + q ln(« p) 



(6.22) 



so that r and p are related by 



p = r 



\ q ln(/tr) 



(6.23) 



Hence, to the order considered, GMj r w GMj p. Using (|4.18J) we obtain the following 
result for the metric in the standard r-coordinate system: 



ds' 



, 2GM , / \ 
1 h q ln(«r) 

r 



dt 2 + 



2GM 
1 + q 



dr 2 + r 2 da 2 . 



(6.24a) 
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The corresponding position dependence of Newton's constant is given by 



G{r) = G[l +jV(r)] (6.25) 

with 

Af(r) = q ln(«r). (l6~2lb ) 

The result (j6.24l a.b) is valid for L, GM <C r <C L exp(l/g). To lowest order, for 
q = 0, eq. (J6.24a|) equals the linearized Schwarzschild metric. The terms oc q are exactly 
those which we found in Appendix [X] by explicitly solving the linearized field equations. 
In fact, ()6.24a|) is precisely the special case of (jA.19J) which follows from the M- function 

Within its domain of validity, the line element ()6.24a|) is close to the Minkowski 
metric. It is the counterpart of (J5.13)) which we had derived for the fixed point running 
G(k) oc 1/k 2 . Like there we can read off a potential from g u = — [l + 20 m N] which 
governs the modified Newtonian limit: 

GM q , . . . 

mN r = + Jln(«r). 6.26 

r 2 

While in the fixed point scenario the extra term in m N was growing oc r 2 , we observe a 
much milder logarithmic growth in the power law model with a small q. 

The time-like geodesies of massive test particles at nonrelativistic velocities can be 
found from (J6.26)) using the Newtonian equation of motion x = — V</> m N- It implies 
■ (r) for the Newtonian velocity v = r (d<p/dt) along circular orbits. For ()6.26J) 



v = r 



'mN 



we obtain 

s GM q 

v 2 (r) = + 1 6.27 

r 2 

displaying the limiting velocity v 2 (r — > oo) — q/2 = v 2 ^. Eq. ()6.27j) coincides with our 
earlier result (|6.15J) since GM/p differs from GM/r only by a term of second order which 
we neglect. 

When the metric (j6.24aj) is valid the curvature of spacetime is small and an es- 
sentially euclidean interpretation applies, with t playing the role of "the" time, and r 
differing only slightly from the radial proper distance. The assumptions which led us 
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to (|6.24|) guarantee that the modified Newtonian potential is always much smaller than 
unity, m N <C 1. One of the consequences is that the proper time r measured by a clock 
at rest, i. e. with constant coordinates (r, 9, (p), differs from the coordinate time t only by 
a factor dr/dt = 1 + m N( r ) which is very close to unity Thus, within the precision of 
our calculation, the velocity v 2 of (|6.27|) coincides with 2nr divided by the period r of 
the orbit as measured by a clock at rest: v = 2nr/T = (27rr/r) [l + O (0 m N)l ■ Note 
that using the standard form of the metric the proper circumference of a circle of radius 
r equals precisely 2irr. Note also that, according to ()6.23j) . r differs from p by a term of 
order m N only. Hence, within the present approximation, we may identify v with 2nr/T, 
or with 27rr/r, or also with 2irp/T. This justifies our use of the Newtonian language when 
we interpreted the rotation curves in Subsection 16.31 

6.5 Galactic scales revisited 

The results of the previous subsections indicate that the scale-free power law G oc 
k~ q could be able to explain at least the gross features of the "non-Keplerian" dynamics 
of galaxies without the need for dark matter. For a point mass it generates a modified 
Newtonian potential with a ln(r)-term, implying a flat rotation curve. 

Contrary to the IR fixed point scenario of Section the k~ q -mode\ trajectory has 
no direct relation to the RG-improved cosmology of [17,18]; in particular, L = k" 1 and q 
are free parameters, so far unconstrained by cosmological observations. 

Nevertheless the general physical picture is still the same as sketched earlier: Be- 
tween laboratory and cosmological scales, G(k) "crosses over" from a G = const-regime 
to the fixed point regime G(k) oc 1/k 2 . In between, at the scale of galaxies, G(k) already 
shows some nontrivial running, it increases with increasing distance, but the precise form 
of its k- or r-dependence we do not know yet. The interpolating ansatz (|5.1|) seems 
problematic, while ()6.1j) works much better. 

For explaining the rotation curves it suffices that the fc _9 -model trajectory (|6.1|) is 
approximately correct near galactic scales only. We may think of it as part of a trajectory 
G{k) oc k m which interpolates between G = const and G(k) oc k~ 2 by "adiabatically" 
shifting the anomalous dimension t?n from n^ = to its value at the fixed point, 7/n = —2. 
All we need for our above derivation to apply is a range of k- values, comprising 2 or 3 orders 
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of magnitude say, for which = —q is approximately constant, implying v(r) ~ const 
between 1 kpc and 100 or 1000 kpc, say. 

In principle the numerical value of q can be calculated by solving the RG equa- 
tion and identifying the specific trajectory which is realized in Nature. (As in QED, 
for instance, this can be done by fixing certain "renormalized" parameters according to 
their experimental values.) If the trajectory Nature selects displays an intermediate k~ q - 
behavior with q of the right order of magnitude, we would predict that the m N of any 
point mass contains a nonclassical term q In r where the prefactor q is universal, i. e. in- 
dependent of the mass M. Hence the predicted plateau velocity = y/q/2 is the same 
for all masses. 

Within a set of typical spiral galaxies the measured v oo-values differ by a factor 
of 2 or 3, say. While the plateau velocity of real galaxies is not perfectly universal it 
does not vary by orders of magnitude on the other hand, so our prediction of a universal 
Voo might be a sensible first approximation. In reality there exists a correlation between 
the maximum rotation velocity f max and the luminosity L of a galaxy, the Tully-Fisher 
relation t>^ ax oc L. With t> max = v^, and since M/L is approximately constant for spirals, 
this means that Voo has a certain mass dependence: Voo cx M _1//4 . 

The term q lnr which we derived from the /c _<? -running is similar to the MOND 
potential (|1.10|) . The only difference is its prefactor oc \fM which replaces our constant 
q and produces the desired oc M -1 / 4 . In fact, the (nonrelativistic!) field equation 
of the MOND approach [31] is designed in precisely such a way that its solution for the 
potential of a point mass contains a logarithm with a prefactor proportional to \f~M . 

Since in many applications MOND is a successful phenomenological theory it is 
plausible to conjecture that the exact Quantum Einstein Gravity dynamically generates 
the 0MOND-term (jl.lUj) or something similar as a quantum effect in the Newtonian regime. 
Clearly a proof of this conjecture is extremely difficult; but on the other hand it is quite 
easy to understand why our above approximate analysis is not precise enough to reveal 
the comparatively subtle M _1//4 -dependence. 

The RG trajectory, computed exactly or in a truncation, is a universal object which 
does not "know" about the concrete system under consideration. In particular the function 
G = G(k) is the same for all galaxies and independent of M therefore. A M-dependence 
can enter only via the cutoff identification k = k M {r) which makes GmW = G{k = 
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&m(?")) a potentially M- dependent function of the distance. For the cutoff k = £/d g this 
is indeed the case in general because g^ u depends on M, but in the Newtonian limit this 
M-dependence drops out and we are left with k oc 1/r, implying a M-independent G(r) 
and In r-term. In the present approach we can hope to find the Tully-Fisher law only once 
we have gained a much better understanding of the physical cutoff process [15] and the 
resulting refined cutoff identification which, most probably, will deviate from k oc 1/r at 
subleading order. 

It always had appeared problematic to construct a relativistic generalization of 
MOND theory [31]. Leaving the Tully-Fisher issue aside, we can say that the improved 
action approach based upon a power law running of G does indeed provide a relativistic 
generalization of this kind. 

We shall return to the discussion of galaxies in the RG framework in Subsection 17. 71 

7 The modified field equation in the Newtonian regime 

In the previous sections we constructed exact vacuum solutions of the modified 
Einstein equation and we found that, in the Newtonian regime already, those solutions 
display features which suggest that the observed rotation curves could be interpreted as 
a renormalization effect of G. For this reason we are now going to analyze the Newtonian 
limit of the modified Einstein equation itself rather than looking at special solutions. As 
a by-product, this discussion will make it clear how the running-G scenario relates to the 
"standard speculation" about the origin of the flat rotation curves: dark matter. 

Contrary to the previous sections we shall now allow for the presence of "ordinary" 
matter, described by an energy-momentum tensor T^ u of the perfect fluid type. It will 
not be traceless in general so that the results of this section cannot be obtained using the 
Weyl technique. 

Except for staticity we shall not impose any symmetry requirements. Apart from 
questions related to the cutoff identification the discussion will also apply to systems 
which are not spherically symmetric. 
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7.1 Background, metric ansatz, small quantities 

Newton's constant is assumed to be time-independent but to have a weak position 
dependence. We use the parameterization 



G(x) =G[l+JV(x)] 



(7.1) 



and assume that A/"(x) <C 1. Being interested in small isolated systems, the cosmological 
constant is assumed to vanish, A(x) = 0. 

As in the standard discussion of the Newtonian limit we make an ansatz for the 
metric of the form 



where it is understood that 0(x) <C 1 and cj(x) <C 1. Here and in the following we use 
cartesian coordinates x M = (t, x l ) = (t,x) and write dx 2 = Sij dx 1 dxK 

When evaluating the field equations we are going to assume that 0, u, and M are 
of the same order of magnitude. We retain terms of first order in 0, u, M <C 1 and their 
derivatives, but discard higher powers (0 2 , • ■ • ) and cross terms (<f) A/", • - ■ )■ 

In the spacetime (|7.2j) . test particles with nonrelativistic velocities v 2 = (dx/dt) 2 
1 obey a Newtonian equation of motion. Retaining only terms linear in 0, c<j,v 2 , the 
geodesic equation boils down to 



Obviously <p is the "modified Newtonian potential" , denoted m N in the previous sections. 
To first order in <p and us the Ricci tensor of the metric ()7.2|) is given by 



Here V 2 = <5 lJ didj denotes the ordinary Laplacian on flat 3-space. Likewise the compo- 
nents of the Einstein tensor are 




(7.2) 



dt 2 



x(t) = -V0(x(t)). 



(7.3) 



R 00 = V 2 0, R 0i = 0, Rij = did,j (u-4>) + S i:j V 2 u. 



(7.4) 



oo — 



(7.5) 
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7.2 



The tensors 6^ and AT^ 



It is a particularly attractive feature of the modified Einstein equation that, in the 
Newtonian limit, it becomes completely independent of the tensor Q^ v whose form was 
not completely fixed by general principles. In fact, every tensor 8^ which is at least 
quadratic in the derivatives d^G = Gd^Af is of order Af 2 and can be neglected therefore. 
For example, 6^ u = e is indeed bilinear in d^G, see eq. (I2.17j) . Since it is impossible to 
satisfy the consistency condition identically with a tensor which has less than two factors 
of d^G we conclude that 9^ v = + higher orders. 

Upon evaluating (J2.7j) for the metric (|7.2j) we obtain for the components of At^ u : 



At 



oo 



-V 2 W, At Ql = 0, At 



did j 



Sy V 2 



Af. 



The corresponding components of AT^ U are 

AT 00 = - (8ttG) _1 V 2 Af, 
AT 0i = 0, 



ATij = - (8ttG) 



didj - 5ij V 2 



Af. 



(7.6) 

(7.7a) 
(7.7b) 
(7.7c) 



The expressions (|7.6j) and ()7.7j) are correct up to terms of order Af 2 , Af '<f>, and Af lv. 



7.3 Content of the consistency condition 

The field equations whose Newtonian limit we are interested in consist of the mod- 
ified Einstein equation with A = and 6^ = 0, 

= 8ttGT^ + AV, (7.8) 

coupled to the consistency condition which is most conveniently analyzed in the version 
of flUnp . In fact > since D ^ G = At^ = 0(Af), and R pv = 0{w,<f>), its first 

term oc D^G (At^ — R pv ) is of higher order and can be neglected. The same is true 
for the second because = 0(Af 2 ), and the third one vanishes since A = 0. What 
remains is the condition 8n GTJ 1 d^Af = O (Af 2 ,Af(f),Afu). Since Af is time-independent 
but an arbitrary function of x this condition requires that T v % = O (Af 2 , Af<fi, Afu). In 
the metric (|7.2|) the energy-momentum tensor of a perfect fluid has the structure TJ 1 = 



48 



diag [— p, p, p, p\- Imposing T v l = + higher orders implies p = + higher orders and leaves 
p unconstrained. 

Thus we conclude that in the Newtonian limit only pressureless matter can be 
coupled consistently. Henceforth we shall employ the energy- momentum tensor of dust, 
TJ 1 = p5®5q. Note, however, that we are not neglecting the pressure terms in AT^. 

7.4 From Einstein's to Poisson's equation 

When we insert ()7.5|) and ()7.6|) into the modified Einstein equation ()7.8|) we see that 
it decomposes into the (OO)-component equation 

V 2 (w + iJV) =4ttG(x)p (7.9a) 

and the (ij)-components 

(cu - + W) = (7.9b) 

where i,j = 1,2,3. The (Oi)-equations are satisfied identically to the relevant order. 

The system f!7.9|) can be solved in a very elegant way. Eq. (j7.9b|) tells us that 
uo — cf) + M must be constant, and since we would like to recover the standard result oj = <f> 
if M = 0, this constant must be zero: u — <f> + J\f = 0. This allows us to express u in 
terms of the (unknown) field and the (prescribed) function A/": 

uj = <p-H. (7.10) 

Inserting (|7.10j) into ()7.9a|) we obtain a modified Poisson equation for (f>: 

V 2 [<f>-\M] =47rG(x)p. (7.11) 

This equation differs from that of classical gravity in two ways: on its LHS the potential 
is replaced with — AA/2, and on the RHS G is allowed to depend on x. 

The equations (|7.1()jl and ()7.11j) show that the metric of the modified Newtonian 
limit can be found in the following way. Given the position dependence of Newton's 
constant, A/"(x), and the matter energy density p(x) one first solves the Poisson equation 

V 2 = 47rG(x)p = 47rG(l+AT(x))p. (7.12) 



didj - 8ij V 2 
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Then, knowing the auxiliary potential 0(x), one computes the metric components as 

(7.13a) 

u = $-\M. (7.13b) 

The resulting line element reads 



ds' 



1 + 26 + Af 



dt 2 + 



2cj) + M 



dx 2 



(7.14) 



Eqs. (|7.12j) and ()7.13a|) are the most important result of this section. They show 
that the modified Newtonian potential is given by 



(x) + ±A/"(x), 



(7.15) 



where is a solution to the "improved Poisson equation" (|7.12j) which involves a position 
dependent Newton constant. Adopting natural boundary conditions, 



d3x , G(xQ P (xQ 



(7.16) 



Eq. ()7.15|) is a generalization of the result (jA.17a|) which we derived in Appendix 
El for the spherically symmetric case using a completely different method. Note that the 
spatial part of the line element (|7.2j) differs from that of the linearized standard metric 
()A.1|) employed in the Appendix. 

Up to terms of higher orders, the metric (|7.14|) admits the factorization 



ds 2 



;i+ao 



dt 2 + 



d^ 



(7.17) 



This factorization is reminiscent of the Weyl transformation (J4.2j) . In fact, if p = for 
|x| > then 0(x) = — GM/\x\ with G = G(x = 0). As a result, the solution belongs 
to Class II, and (j7.17j) is the familiar Weyl rescaling ds 2 = [G(x)/ G] ds 2 applied to the 
linearized Schwarzschild metric. For p ^ the result ()7.17|) with (j7.16J) is new and cannot 
be proven by means of the Weyl trick since T ^ ^ in this case. 

For a 5-function source at the origin equals the classical 1/r-potential, and ac- 
cording to ()7.15|) the corresponding modified Newtonian potential is 



(x) = -^ + -Af(x). 



(7.18) 
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The generality and simplicity of this result and its generalization ()7.15|) is striking: 
The impact of the x-dependent G on the Newtonian limit consists of the additional term 
A/"(x)/2 in the potential only. This is true for any externally prescribed weakly position 
dependent Newton constant. The function G(x) is neither required to have any special 
symmetry nor is it necessary that G(x) derives from a running G(k) via some cutoff 
identification k = k(x). 

While this is difficult in general [15], it is fairly straightforward to convert G{k) to 
a position dependent G if the system is spherically symmetric. In the Newtonian limit 
this step simplifies considerably. In the relevant order the cutoff identification k = k(r) 
is simply k = £/r, r = |x|, since when inserted into Af(k = £/r) = Af(r) the difference 
between r and d g (r) amounts to a higher order term. 

As a result, the potential of a point mass at the origin is given by <p(r) = —G M/r + 
M{r)/2. The renormalization effects generate an additional contribution —N'(r)/2 to the 
classical force —GM/r 2 . The new force is attractive if jV(r), and hence G(r), grows as a 
function of r. 

In Sections E] and IH1 we encountered two examples already. There we calculated the 
exact metrics pertaining to the fixed point and the power law running of G(k), respectively, 
and we saw that in certain parts of those spacetimes a Newtonian interpretation is possible. 
By eqs. f|5.13|) and ()6.24b|) the corresponding A/"-functions are 

jV(r) = k 2 r 2 (fixed point) (7.19a) 

jV(r) = q ln(/«r) (power law) (7.19b) 

Using these TV's in eq. (|7.18|) we reproduce exactly the potentials (|5.14|) and (|6.26|) ex- 
tracted from the exact solutions. In Sections El and El we had to assume that 9^ u = O^P. 
Now we see that the modified potentials obtained there are actually independent of the 
choice for the ^-tensor. 

7.5 The importance of the vacuum pressure 

Using the explicit expression for AT 00 , eq. (|7.7a|) . we can rewrite the Poisson equa- 
tion (|7.11|) in the following equivalent way: 

V 2 = 4ttG(x) [T o-AToo]. (7.20) 
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At first sight it seems quite puzzling why the negative of AToo, the "vacuum" energy 
density due to the x-dependence of G, should be added to the matter energy density 
p = Too- This paradox is resolved if one recalls that even in the Newtonian limit of 
standard general relativity [43], when the pressure is not negligible, Newton's T Q0 gets 
replaced by T 00 + Ylt=i Tu (which equals the familiar p + 3p for a perfect fluid): V 2 = 

47rG[T 00 + £ti T J- 

Above we neglected the pressure of the ordinary matter but the pressure terms 

stemming from AT M „, i.e. AT^-, where fully taken into account. The AT^-components 

f!7.7jl obey the following identity which is quite remarkable: 

3 

AT 00 + ^AT 4i = -AT 00 . (7.21) 

i=l 

Adding the pressure terms has precisely the effect of flipping the sign of AT 00 . Therefore 
we may rewrite (|7.2(J|) in the form 

3 

V 2 = 4tt G(x) [Too + AT 00 + ^ AT,,] . (7.22) 

i=\ 

The source term of this Poisson equation has exactly the expected "p + 3 p" -structure and 
we understand that there is nothing wrong with the difference Too — AToo appearing in 
(I7.20J) . In fact, this relative minus sign resulting from the inclusion of the vacuum pressure 
is crucial in order to obtain a physically sensible Newton limit. If we naively use the sum 
of the energy densities T o + AT o as the source for </>, the ensuing modified Newtonian 
potential is given by (j7.18J) with A/"(x) replaced by — JV(x). The consequence would be 
that a G(r) growing with r gives rise to a repulsive component in the gravitational force. 
This is inconsistent with the antiscreening nature of gravity and indeed is not what the 
theory predicts. 

7.6 Effective energy density vs. dark matter 

For the interpretation of the RG improved theory it is helpful to rewrite the modified 
Poisson equation in the style of a classical Poisson equation involving a modified energy 
density, though. Defining the effective density 

p eff = [G(x)/G] (p-Ap), (7.23) 



52 



eq. ()7.20|) assumes the traditional form 



V 2 = 47rGp eff . 



(7.24a) 



Above we introduced the vacuum energy density 



Ap = AT C 



oo — — 



(7.25) 



in analogy with T 00 = p (to lowest order). Since Ap is of order Af already, (|7.23|) is 
equivalent to 



What is the physical significance of p e ff? 

In astronomy a standard method for measuring mass distributions proceeds, schemat- 
ically, as follows. One monitors the trajectories x(t) of as many "test particles" as possible, 
assumes they satisfy x = — V0(x), then reconstructs from the trajectories, and finally 
computes the matter density responsible for this potential as (V 2 0) /47rG, taking the 
Newtonian field equation for granted. 

In the Newtonian regime, the RG improved theory of gravity proposed in this paper 
consists of two parts: the (modified) field equation (|7.24|) and the (standard) particle 
equation of motion ()7.3|) . Let us assume both of these laws are correct, and that they do 
indeed govern the dynamics of "test particles" in a galaxy, say. Then, when an astronomer 
tries to measure the mass density inside the galaxy by applying the method described 
above, he or she will come to the conclusion that this mass density is given by p e fj rather 
than the density p of ordinary matter. In particular, in a halo with p = the measurement 
would determine p c s = —Ap. 

The most popular speculation about the origin of the flat galaxy rotation curves 
is the assumption that there exists a certain amount of dark matter in galaxies, with 
density Pdm> which generates an appropriate potential in the conventional way, i.e. via 
V 2 = 4:7iG (p + Pdm)- Comparing this equation to ()7.24|) we see that in the RG approach 
it is essentially — Ap which takes the place of Pdm- in fact, typically the AAp-term in 
()7.24bj) is tiny compared to p so that p c g « p — Ap. 

Thus we conclude that a scale- and, as a consequence, x-dependent Newton constant 
can mimic the presence of dark matter. In the Newtonian limit the effective density of 
the dark matter substitute is given by +(8ttG ) V 2 G(x). 



p eS = P + Mp+ (8ttG) V 



2 M. 



(EM)) 
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7.7 A simple model of a galaxy 

Let us consider a spherically symmetric, static system so that p and M depend on 
r = |x| only. Then the effective density is given by 

p cff (r) = [1 +JV(r)] p(r) + (svrGr 2 )" 1 A (V 2 . (7.26) 

The Poisson equation (|7.24a|) is solved by integrating 



A, w = ^ (7 . 27) 



with A1 e fr(r) = An J dr' r' p e s(r'), or explicitly, 
o 

2 j 

^cfr(r) = M(r) + — JV(r). (7.28) 
2G dr 

r 

Here .M(r) = An J dr' r' 2 (l + A/"(r')) p(r') is essentially the mass of the ordinary matter 
o 

contained in a ball of radius r, centered at the origin. Since M <C 1, the factor (1 + A/") 
is very close to unity and plays no role at the qualitative level; we shall ignore it in the 
following. The effective mass inside the ball, A4 e s, consists of the classical part A4, plus a 
term oc djV/ dr which is entirely due to the renormalization effects. If Newton's constant is 
a growing function of r then the nonclassical contribution is positive and -M c fr(?") > -M(r). 

In the jargon of quantum field theory one would refer to M. and Ai e s as the "bare" 
and the "dressed" or "renormalized" mass, respectively. The antiscreening character of 
Quantum Einstein Gravity [5] amounts to cUV /dr > so that the quantum-gravitational 
"dressing" of a classical mass distribution increases the total mass. 

From ()7.27|) with (|7.28J) we obtain the potential 

} GM(r') 1 
0(r) = JdT'^^ + -U{r). (7.29) 

The corresponding rotation curve v 2 (r) = r<ft'(r) for particles orbiting in this potential is 

9/ s GM(r) 1 d ... s ,„ „ . 

v 2 r = — + -r—J\f(r). 7.30 

r 2 dr 

On the basis of the potential (J7.29j) we can make a simple model of a spherically 
symmetric galaxy. We identify p with the density of the ordinary luminous matter. We 
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model the luminous core of the galaxy by a ball of radius tq. The mass of the ordinary 
matter contained in the core is Ai(r ) = Aio, the "bare" total mass of the galaxy. Since, 
by assumption, p = and hence A4(r) = Aio for r > r , the potential outside the core is 

' 1 C M ' + \N{r) (r>r ) (7.31) 



rw r 2 

We shall refer to the region r > r as the "halo" of the model galaxy. For r sufficiently 
large, the potential in the halo is essentially M{r)/2. 

Which fc-dependence of G(k) is phenomenologically favored? We saw that a power 
law behavior G(k) oc k~ q is, on the one hand, the simplest and most natural option 
theoretically since it arises from a constant anomalous dimension. On the other hand 
it was precisely this fc-dependence which led to flat rotation curves for point masses. 
Employing the identification k oc 1/r the r-dependence of G is governed by a power law 
as well: G(r) oc r q . Only if q is small and r is not too large (r <C L exp(l/g)) we may 
expand r q , obtaining G{r) = G [l + jV(r)] with M{r) = q ln(fi;r). 

Accepting the logarithmic A/"-f unction as the most promising candidate, the position 
dependence of G generates the negative vacuum energy density 

87rG r z 

It also generates vacuum stresses (pressure terms) thanks to which the effective energy 
density receives a positive vacuum contribution: 

Peff(r) = p(r) + (7.33) 
87rG t 1 

This formula is valid for all r > r tr where r tr is the radius at which the nontrivial RG 
running gets "switched on" . 

For simplicity we assume that r tr < r , i. e. that the transition takes place inside the 
core. As we discussed already, the scale r tr is comparatively unimportant: the renormal- 
ization effects become strong, and visible, at the larger scale r nc = 2GAio/q only. Thus, 
everywhere in the halo, 

0( r ) = _9^k + | i n ( K r ). (7.34) 

Obviously (J7.34j) has the same form as the potential (J6.26)) for a point mass at the origin. 
Strictly speaking ()7.34j) is valid for ro < r ^ L exp(l/g) only. However, most probably 
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the upper bound on r is irrelevant from the practical point of view since at too large 
distances the picture of a single isolated galaxy is unrealistic anyhow. At a certain point 
the influence of other galaxies and of the expansion of the Universe should be included. 

Of course the most important experimental fact supporting the power law running 
of G are the observed flat galaxy rotation curves. The power law scenario leads to the 
prediction 

2 GM(r) q 
v ( r ) = " + 2 ( 7 ' 35 ) 

which is valid for r > r tr , both inside and outside the core. In the halo we have v 2 (r) = 
GMo/r + q/2 and the situation is the same as for the point mass in Section |H1 For 
r — ► oo, ()7.35|) yields a constant velocity = q/2. As we discussed already, the observed 
values of Voo suggest that q = 0(lCr 6 ). The smallness of this number justifies the 
linearization with respect to N '. It also implies that the variation of G inside a galaxy 
is extremely small. The relative variation of Newton's constant from some r\ to r 2 > r% 
is AG/G = q\\i{r 2 /ri). As the radial extension of a halo comprises only 2 or 3 orders 
of magnitude (from 1 kpc to 100 kpc, say) the variation between the inner and the outer 
boundary of the halo is of the order AG/G ~ q, i.e. Newton's constant changes by one 
part in a million only. 

The r-dependence of the velocity (|7.35j) is in qualitative agreement with the obser- 
vations. For realistic density profiles, Ai(r)/r is an increasing function for r < ro, and 
it decays as Aio/r for r > r . As a result, v 2 (r) rises steeply at small r, then levels off, 
goes through a maximum at the boundary of the core, and finally approaches the plateau 
from above: v 2 (r > ro) = GM.$/r + q/2. Some galaxies indeed show a maximum after the 
steep rise, but typically it is not very pronounced, or is not visible at all. The prediction 
for the characteristic scale where the plateau starts is r nc = 2GM.§/q. With q = 10~ 6 and 
M.q = 10 11 Mq one obtains r nc = 9 kpc which is just the right order of magnitude. 



We close this section with a set of further comments. 

(a) As for deriving more accurate, quantitatively correct rotation curves the limiting 
factor of our method is the cutoff identification k = k(x). Beyond spherical symmetry, for 
a disk say, there will be more than one relevant length scale and there is no straightforward 
way of finding out which function of those scales will take the place of the essentially 
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unique single-scale ansatz k oc 1/r. In Subsection 16.51 we explained already that for a 
very similar reason the Tully-Fisher relation cannot be obtained easily. In principle it is 
clear how to solve these problems: One computes T = T^o and then solves the effective 
field equation 6T/5g^ v = 0. This is a formidable task, though. Only under very special 
circumstances, strong symmetry constraints for example, the improvement strategy is 
available as a technically feasible "short-cut" to the full-fledged calculation. (See [15] for 
a detailed discussion of this point.) 

(b) The actual mass of the model galaxy, including the renormalization effects, is 
M. e s- Integrating (|7.33j) we find for this "dressed" mass 

M eS (r)=M +(q/2G) r (r > r ) (7.36) 

The first contribution, the "bare" mass Aio, is the total mass of the ordinary luminous 
matter contained in the galaxy, the second is due to the running of G. The latter con- 
tribution increases with r unboundedly so that, strictly speaking, the total dressed mass 
A^efr(oo) is infinite. However, since it is clear that our assumptions break down at some 
large value of r the growth of M. c r will stop at some point probably, but this is outside 
the scope of the present model. 

(c) In the dark matter scenario it is usually assumed that galaxies form in two 
steps. In the first step, density perturbations which have grown sufficiently lead to the 
formation of virialized, gravitationally bound structures consisting of dark matter only. 
Then ordinary (baryonic) matter will cool by radiating energy, sink to the centers of those 
dark matter haloes, and thus form galaxies [47]. 

In the physical picture we advocate in the present paper quantum fluctuations of the 
gravitational field take the place of dark matter. If we try to keep the essential features 
of the above structure formation scenario unaltered it is natural to ask which objects 
could play the role of the pure dark matter haloes which contain no ordinary matter 
yet. Among the exact solutions which we found above there is an obvious candidate: 
the A4q = M = special case of the metric ds 2 = w(p) ds 2 . It is generated by the 
Minkowski metric 7^ = 77^, whence g^ u = w(p) r]^ v with the Weyl factor given by ()6.7|) 
with If(p) = p. This metric describes a "solitonic" excitation of the gravitational field 
in absence of any matter. Contrary to the Schwarzschild spacetime which is a vacuum 
solution, too, it owes its existence entirely to the renormalization effects: for q — ► it 
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approaches flat space. 

Note that if one interprets the 4-momentum carried by the vacuum fluctuations as 
due to some kind of matter, this matter is indeed "dark" in the sense that it is not "seen" 
by the electromagnetic interaction, the reason being that the solution g^ v is conformally 
related to rj^. 

It seems conceivable that in the early Universe density fluctuations first generate the 
"solitons" with M — 0, thereby creating gravitational sinks which get filled with ordinary 
matter later on. It would be interesting to see if the current dark matter-based scenarios 
of structure formation can be reinterpreted in a running-G framework. This is beyond 
the scope of the present paper, however. 

8 Summary and conclusion 

As in any quantum field theory also in Quantum Einstein Gravity there is a standard 
way of investigating quantum effects both at small and large distances: one computes 
the effective action V [g^ u ] and then solves the effective field equation 5T/5g^ u = for 
g^v, the expectation value of the microscopic metric. A sufficiently precise approximate, 
let alone exact computation of T [g^ v ] is an extremely difficult task, well beyond our 
present technical abilities. However, if nonperturbative Quantum Einstein Gravity is 
asymptotically safe, and there is growing evidence that this is actually the case, then the 
problem is indeed more of a computational rather than conceptual nature; many of the 
problems encountered in the traditional approaches (perturbative nonrenormalizability 
etc.) would be overcome then [9]. 

Under certain particularly favorable physical conditions there exists a "short-cut" on 
the way to the expectation value g^ v {x) which renders the calculation of Y unnecessary: 
renormalization group improvement. The idea is to use a truncated form of the flow 
equation governing the scale dependence of the effective average action Y k [g^ v ] in order 
to obtain running gravitational couplings G(k), A(k), ■ ■ • which then, in a second step, are 
converted to scalar functions on spacetime by means of a "cutoff identification" k = k(x). 
The most subtle issue of this approach is finding such an identification which requires 
a detailed understanding of the relevant physical scales involved [15]. In ref. [15] we 
proposed to use the scalars G(x) and A(x) thus constructed for an improvement of the 
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Einstein-Hilbert action (as opposed to its field equation or solutions thereof). 

While in ref. [15] we focused on general aspects of the resulting Brans-Dicke-type 
effective theory, as well as on applications to cosmology, the present paper is dedicated 
to spherically symmetric systems. Our investigation of static isotropic solutions to the 
effective field equation in Sections H3 and H] has a variety of applications, only some of 
which we already touched upon. In particular the methods developed there will be the 
basis for a discussion of black holes in the improved action framework, a topic to which 
we shall come back elsewhere. 

In the present paper we analyzed how gravity gets modified at very large distances 
if one assumes that Newton's constant has a nontrivial RG running in the infrared, i. e. 
at astrophysical scales. As far as the pertinent RG trajectories are concerned this issue is 
more speculative than the black hole problem but it can be analyzed more easily since the 
essential physical effects are of a quasi Newtonian nature. As an application we explored 
the possibility that the flatness of the observed galaxy rotation curves, usually attributed 
to the presence of dark matter, is actually due to a RG running of G. 

From the point of view of the present analysis, the "RG trajectory" G(k) is a given 
function to be confirmed - or falsified - by future flow equation studies. The "IR fixed point 
scenario" which we analyzed in Section|5]was inspired by recent work [17,18] which showed 
that the fc-dependence G(k) oc 1/k 2 , at cosmological scales, leads to a phenomenologically 
viable solution of the "cosmic coincidence problem" ; it explains the approximate equality 
of the present matter and vacuum energy density without a finetuning of parameters. We 
employed the model trajectory G(k) = G + g* R /k 2 which interpolates linearly between 
classical gravity in the laboratory and the l//c 2 -running at cosmological distances. We 
constructed exact vacuum solutions for this ^-dependence, analyzed their properties in 
various regimes, and found that the modified Newtonian regime (realized at intermediate 
distances) could perhaps be relevant to the physics of galaxies. Here Newton's 1 / r 2 -force 
is modified by a nonclassical contribution linearly increasing with r. However, it turned 
out that this force is too weak to explain the observations. 

In Section El we maintained the basic idea that Nature interpolates between G = 
const in the laboratory and (something similar to) G(k) oc 1/k 2 in cosmology but we 
employed a different interpolation. One could imagine that the anomalous dimension 
7? N = kdkG/G interpolates smoothly between r/ N = and its fixed point value —2 with 
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a very small rate of change, at least within some interval of fc-values. In particular we 
assumed, and this is the only assumption relevant to the present study, that r] N = —q 
is approximately constant for k ranging over the 2 or 3 orders of magnitude which have 
an impact on the structure of galaxies (corresponding to distances between 1 kpc and 
1 Mpc, say). In this domain Newton's constant scales as G(k) oc k~ q with a small positive 
exponent q. 

Also for the fc _l3 -trajectory we were able to find exact vacuum solutions generalizing 
the Schwarzschild metric, and we analyzed their geodesies and the corresponding rotation 
curves. We found that at intermediate distances there exists a quasi Newtonian regime 
where Newton's 1/r-potential is supplemented by a term q lnr which is well known to 
produce a constant v(r). While this regime might indeed be relevant to the physics 
of galaxies, the full nonlinear vacuum solutions are interesting also in their own right. 
They are based upon a background independent, selfconsistent cutoff identification k oc 
l/d g (r) where, contrary to earlier approaches, the proper distance d g {r) is computed from 
the solution g^ v rather than from a fixed background metric; as a result, the position 
dependence G = G{r) originating from a given G = G(k) is computed rather than put 
in by hand. A particularly remarkable vacuum solution obtains by applying the solution- 
generating Weyl transformation to Minkowski space. It represents a kind of gravitational 
"soliton" which can exist only thanks to the running of G. 

The "phenomenologically" most promising portion of the spacetimes constructed 
turned out to be in the weak field regime. This motivated us, in Section to investigate 
the Newtonian limit of the modified Einstein equation in general, rather than of special 
solutions, with matter included and without insisting on spherical symmetry. One of 
the key results was that any arbitrary G(r) which increases with distance always gives 
rise to an attractive nonclassical force which adds to Newton's l/r 2 -contribution. If 
G — G [l + A/"] , a term +J\f/2 gets added to the classical potential. 

The upshot of all these investigations is that it is very well conceivable that the 
approximately flat galaxy rotation curves one observes in Nature are actually a renormal- 
ization effect of Newton's constant rather than due to the presence of hitherto undetected 
dark matter. If so, the fc-dependence of G(k) in the range of typical galaxy scales is 
governed by a power law most probably. In Subsections 16.51 and 17. 71 we described in detail 
its consequences for the physics in galactic haloes and the origin of the constant rotation 
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curves in particular. We also discussed various limitations of the present calculational 
scheme; for instance, it will be difficult to go beyond spherical symmetry or to describe 
higher order effects such as the Tully-Fisher relation. 

Now that we know which type of RG trajectory we are looking for it will be exciting 
to see whether future RG studies, employing improved truncations, can confirm this 
picture. In a sense, it would replace dark matter with the vacuum fluctuations which 
drive the RG evolution of the gravitational parameters. 



Appendix 

A Linearized static isotropic field equations 

In this appendix we analyze the coupled system of equations (|3.13|) in the modified 
Newtonian limit [43]. For the metric we make the ansatz 

ds 2 = - [1 + 2 0(r)] dt 2 +[1-2 X {r)] dr 2 + r 2 da 2 , (A.l) 

which amounts to 

A(r) = l-2x(r), B(r) = 1 + 2 4>(r). (A.2) 

Deviations from flat Minkowski space are assumed tiny, and the functions <fi and x are 
considered small compared to unity therefore. In the following calculations we retain 
terms linear in 0, x, an d their derivatives, but discard all higher powers or cross products 
of them. As a result, Einstein's equation ()3.10|) consists of the linearized tt-component 

4 + - = -47rGp tot , (A.3a) 

the rr-component 

^- + ^- = 47cGp t r °\ (A.3b) 

and the 00/</?</?-cornponent 

- + - + 0" = 8nGpf\ (A.3c) 
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The consistency condition assumes the form 

£ Pl ot + 0' (p tot + p") + \ (vT - pT) = o, 

and the linearized form of the ordinary continuity equation D^T^ = reads 

p + (p + p) (j)' = 0. 
Next we must linearize the components of AT^. One finds 
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Inserting the above leading order approximations into Einstein's equation (jA.3|) we obtain 
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For 4> and x to be small it is necessary that the sources of the gravitational field 
are weak. Likewise we are going to assume that the variability of G(r) is small since it is 
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obvious from the above equations that also the position dependence of G acts as a source. 
We use the parameterization 

G{r) = G{l+N{r)] (A.10) 

and assume that M <C 1 is of the same order as <fi and \. We retain terms linear in 0, 
X, A/", and their derivatives, but neglect higher powers and products. Hence, to leading 

order, G'/G = M' + • • • , {G'/G) 2 = + • • • , G"/G = N" H , and G'G"/G 2 = + • • • . 

Furthermore, we neglect the cosmological constant, A(r) = 0, because it is not important 
in the present context. As a result, the tt-component (|A.8aJ) boils down to 

* + * = _ 47r G r p + i_ + ^_ (A.lla) 
and the rr-component (JA.8bj) becomes 

Y </>' A"' 

4 + - = 4vrG(r)j9 + — . (A. lib) 

Rather than the 99 / ^(^-component, we consider the consistency condition along with 
the continuity equation the independent equation. Linearizing the consistency condition 
(|A.9jl yields the statement that, up to higher orders, 8n G Af'p = 0. Since, by assumption, 
M' 7^ 0, we obtain the condition p = 0, up to higher orders. This shows that metrics of 
the type (|A.1|) can describe the modified Newtonian limit at most for pressureless matter 
("dust"). 

The last equation we must satisfy is the continuity equation (jA.5|) . It does not con- 
tain M i and coincides with the one appearing in the standard Newtonian limit therefore. 
For p = it says that p<p' = up to terms oc 4> 2 ,x 2 ,(pX- As in the standard case, we 
consider the product pep' negligible, formally of order 2 , so that (jA.5j) is satisfied with 
sufficient accuracy. This is the Newtonian approximation of "quasistatic dust" [45]. 

It remains to solve (jA.lla|) coupled to (jA.llbj) with p = 0. The most elegant method 
is to make an ansatz 

0(r) = 0(r) + ijV(r) (A.12a) 
Y(r) = x(r) + \ r N\r) (A. 12b) 
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with two new functions, and x, yet to be determined. Inserting this ansatz into (|A.lla|) 
and (jA.llb|) for p = yields, respectively, 

^ + ^ = -47rG(r)p, (A.13a) 

4 + - = 0. (A. 13b) 

The equations (|A.13l a.b) have a strikingly simple interpretation. They are nothing but 
the tt-and the rr-component, respectively, of the field equation 

G, v (g)=MG{r)T llv (A.14) 

evaluated for the following metric analogous to (jA.lj) : 

ds 2 =g^dx IM dx u = -[1 + 2 $(r)] dt 2 + [l -2x0)] dr 2 + r 2 da 2 , (A. 15) 

Obviously (jA.14j) is precisely the Einstein equation the improved equation approach is 
based upon: it obtains from the classical = SirGT^ by the substitution G — > G(r). 
While and x determine the metric in the improved action setting, and x P^y the 
analogous role when the field equation is improved. 

The system (|A.13l a.b) is easily solved. (jA.13b|) tells us that x is known once we know 
(J): x = —r (f)' . Inserting this into (jA.13aj) yields the one-dimensional Poisson equation 

VL0 = 4vr G(r)p(r), (A.16) 

where V^ ad = '§^ + ^^ = }:'i^ r ^ s the radial part of the Laplacian on flat space. 

Let us summarize. If the matter is quasistatic dust of density p, the modified 
Newtonian limit implied by the improved action functional is given by the metric (jA.l|) 
with 

<f>(r) = 0(r) + |AT(r) (A.17a) 
x (r) = -r 0'(r) + \ r M\r) (A. 17b) 

where M{r) parameterizes the variation of Newton's constant via (|A.10p . and <p is to be 
obtained by solving the "improved Poisson equation" (jA.16|) . Note that the final result 
does not depend on e, i. e. on the 6*-tensor chosen. 
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The possibility of generating solutions of the modified field equations from solutions 
of a much simpler equation, here the improved Poisson equation, is reminiscent of the 
Weyl transformation technique. However, the relations (jA.17|) just discovered apply to 
matter with = diag[— p, 0, 0, 0]. Hence T = — p ^ 0, so that the solutions are not in 
Class II or Ilia where the Weyl method is available. 

It is instructive to study the special case of vacuum solutions in the classical limit. 
The solution to Poisson's equation, for r ^ 0, is Newton's potential then: 

?(r) = -— = Mr) (A.18) 
r 

In this case — r 0' = = 0n so that x = 0n + \ tM' and the metric assumes the form 

ds 2 = - [l + 2 N (r) + Af(r)] dt 2 + [l - 2 </> N (r) - rU\r)] dr 2 + r 2 da 2 . (A. 19) 

The special case (|A.19|) does indeed belong to Class II, and it can be obtained by a 
Weyl transformation. In the Newtonian limit can be read off straightforwardly from 
ds 2 = [G(p)f G~\ ds 2 because the difference between a selfconsistent cutoff based upon 
d g (p) and the naive one, employing the Minkowskian d-y(p) = p, is of higher order and 
negligible therefore. Hence G(p) = G(k = £//?) = G [l +jV(p)] is independent of g^ u so 
that the explicit solution reads, in the p-coordinate system, 

ds 2 = [1 + H{p)] (-f(p) dt 2 + f(p)- 1 dp 2 + p 2 da 2 ) (A.20) 

with f(p) = 1 + 2 0n(p). Retaining only terms linear in 0n and Af, this line element can 
easily be cast into the AB-iorm. Upon defining the new radial coordinate 

r = p^TTW = p[l + \N +0(N 2 )] (A.21) 

one finds precisely eq. (|A.19|) as the final result. It amounts to the solution 

. , 2 G Ad . . . , , . 

A(r) = l + rj\f (r), A.22a 

r 

B(r) = 1 - + Mr), (A.22b) 

r 

G(r) = G[l+M(r)], (A.22c) 
which is valid for an arbitrary (but small) function M{r) . 
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B Geodesies 



In this appendix we study the trajectories of test particles in spacetimes with a 
metric of the type 



ds 2 = w(p) 



-f(p)dt 2 + f(p)- 1 dp 2 + p 2 (d9 2 + sin 2 9 dp 



(B.l) 



where w(p) and f(p) are arbitrary functions for the time being. Geodesies x M (A), parametrized 
by the parameter A, can be obtained from the Lagrangian 

L = -w{p) f(p) i 2 + w{p) f(p)- 1 p 2 + w(p) p 2 (O 2 + sin 2 9 2 ) . (B.2) 

In this appendix, a dot denotes the derivative with respect to A, and a prime a derivative 
with respect to p. The Euler-Lagrange equations implied by (jB.2|) for t, if, p, and 9, 
respectively, read 

w ft = const = a, (B.3a) 
w p 2 sin 2 9 (p = const = (3, (B.3b) 

2jp+(^ p 2 = -{wf)'i 2 +(w p 2 ) ' (9 2 + sin 2 9 2 ) , (B.3c) 

wp 2 9+ (wp 2 J 9 2 = wp 2 sin9 cos9tp 2 . (B.3d) 

Without loss of generality we may assume that the motion takes place in the plane 
9 = ir/2. Eq. (jB.3d|) is satisfied identically by the ansatz 9 (A) = tt/2 = const, and the 
remaining equations are 

t = a/(wf), (B.4a) 
= /3/( W p 2 ), (B.4b) 

2 j P + (j ) ' P 2 = - (w /)' i 2 + (w p 2 ) ' v 2 . (B.4c) 



We are not going to embark on a general analysis of the equations flB.4|) here. Instead 
we focus on the question whether or not the spacetime admits, for massive test particles, 
closed circular orbits at p = const. If they exist, w f and w p 2 are constant for a given 



67 



orbit, and the equations (jB.4a|) and (jB.4b|) simply tell us that t(A) oc A and y?(A) oc A. 
Inserting (jB.4a|) and (|B.4b|) into ()B.4c|) with p = p = we see that circular orbits do exist 
if constants a and (3 can be found such that 



As expected, there are no solutions in Minkowski space (w — f — 1), but there exist 
circular orbits in the Schwarzschild geometry (w = 1, f(p) = 1 — 2GM/p), say. If one 
specializes for the setting of Section f(p) is fixed to be of the form ()4.3b|) so that it 
is the Weyl factor w(p), related to the p-dependence of G, which then decides about the 
existence of circular orbits. 

The angular velocity referring to A, 6, as well as the angular velocity defined in 
terms of the coordinate time t, dip/dt = <p/t, are constant for a given circular orbit, but 
in general they will depend on p, of course. From (|B.4c|) we obtain 



This quantity is indeed positive precisely if (|B.5|) is satisfied. It determines the angular 
velocity for which the centrifugal force balances the gravitational attraction. The period 
of an orbit at coordinate radius p, measured in units of the coordinate time t, is given by 
T(p) = 2tt/ (dip/dt), or 




(B.5) 




w'f + wf 
w' p 2 + 2 pw 



(B.6) 



2 , «/ 
p w 



,i -i 1/2 



T(p) = 27ip 




(B.7) 



In particular for f(p) given by (Obi . 



1 + Y(p) 



I 1/2 



T(p) = 27rp 



(GM/p)+f( P )Y(p) 



(B.8) 



with the abbreviation 



Y(P) 



p w'(p) 
2 w(p) ' 



(B.9) 



In the classical case, Y = 0, and we recover Kepler's law T 2 
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